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ABSTRACT 


The effect of changes in input variables on the calculated 
bottom-hole pressures of flowing gas wells for six different calcu- 
lation methods were investigated. The six methods studied were: 

1. Average Temperature and Average Compressibility 

Factor Method 

2. Average Density Method 

3. Cullender and Smith Method 

4. Poettmann Method 

5. Sukkar and Cornell Method 

6. Dranchuk and McFarland Method 
Both the analytical and numerical approaches were used. Numerical 
solutions and computer programs were developed for all the six methods. 
However, analytical solutions were developed only for methods 1 and 2. 

The input variables studied were well depth, flow rate, well- 
head pressure, bottom-hole temperature, wellhead temperature, gas 
gravity and absolute pipe roughness. The error introduced by using 
linear interpolation of tabulated integral values in the Poettmann and 
in the Sukkar-Cornell Method was also studied. 

For methods 1 and 2, the results obtained using the analytical 
and numerical approach are in good agreement. For all the methods 
studied, the results showed that the calculated bottom-hole pressure 
is most sensitive to variations in wellhead pressure and least sensi- 


tive to temperature variations. The error introduced by using linear 
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interpolation in the Poettmann Method in the pressure and tempera- 
ture range studied in this work is small and thus linear interpola- 
tion can be used with confidence when using this method. However, 
this procedure is not recommended for use in the Sukkar-Cornel] 


Method, for the error introduced is significant. 
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CHAPTER I 


INTRODUCTION 


Bottom-hole pressure can be determined either by measurement 
with a bottom-hole pressure gauge or by calculation. Since measurement 
with a bottom-hole pressure gauge is costly, especially for deep wells, 
calculation of bottom-hole pressure is to be preferred when possible. 
Most of the calculation methods presently available are based on the 
mechanical energy balance or the general flow equation. Varying the 
assumptions in solving the flow equation has led to many different 
bottom-hole pressure calculation procedures. In all the methods, well 
effluent composition or gravity, wellhead pressure, flow rate, well 
depth and well bore temperatures are factors that must be known before 
calculating the bottom-hole pressure. 

The most widely used and recently developed methods are: 

1. Average Temperature and Average Compressibility 

Factor Method 

2. Average Density Method 

3. Cullender and Smith Method 

4, Poettmann Method 

5. Sukkar and Cornell Method 

6. Dranchuk and McFarland Method 

For Methods 1 and 2, a closed-form solution is available and 
thus analytical methods were developed to solve the flow equation for 
these two cases. For the remaining four methods, Leibniz's formula 


may be used to differentiate under the integral sign. If this 
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aiiiaeeh is taken, analytical equations can be developed for these 
four cases as well. However, the resulting equations will still ae 
to be integrated and, consequently very little is gained by taking 
this approach for methods 3 to 6. Therefore numerical solutions and 
computer programs were developed for all the six methods to calcu- 
late the bottom-hole pressure. 

All the six methods mentioned were originally developed for 
use up to a reduced pressure of 15 or less due to the limit of 
availability of the compressibility factor. This study has extended 
all the above methods to calculate bottom-hole pressure up to a 
reduced pressure of 30, using a recently developed, extended compress- 
ibility factor calculation technique. Gases containing contaminating 
HoS and CO» were handled by applying the Wichert and Aziz [1] pseudo- 
critical properties modification method. 

Very often the bottom-hole pressure calculated using the above 
methods differs from the measured value. This may be due to (a) the 
shortcomings of the method itself, and (b) inaccuracy in the input 
data required to do the calculations. The objective of this investi- 
gation is to examine the effect of (b) for all the six methods. 

The input variables studied were depth, flow rate, wellhead 
pressure, wellhead temperature, bottom-hole temperature, gas gravity 
and absolute pipe roughness. The effect on the calculated bottom- 
hole pressure of varying each variable was studied for all the methods 
mentioned. The effect of absolute pipe roughness was studied for 
Methods 1, 2, 3 and 4. A range of pipe roughness was used as suggested 


by Cullender and Binkley [2], Smith and Williams [3] and Smith and 
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Miller [4]. The friction factor correlation used in Methods 5 and 6, 
which assumes a constant pipe roughness of .0006 inch, is a function 
of pipe diameter only. Consequently, no attempt has been made to 
Study the effect of pipe roughness for these two methods. The effect 
of using linear interpolation in the Poettmann and Sukkar-Cornel] 


Method was also examined. 
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CHAPTER IT 


LITERATURE REVIEW 


A variety of methods has been developed for calculating flowing 
bottom-hole pressure in gas wells. 

R.V. Smith [5] derived an equation for the vertical flow of gas 
Similar to the Weymouth equation [6] for horizontal flow by assuming a 
constant effective compressibility factor and constant average temper- 
ature. This method is known as the Average Temperature and Average 
Compressibility Factor Method. 

Rzasa and Katz [7] assumed a constant average pressure, temper- 
ature and compressibility factor in calculating static pressure 
gradient in gas wells. This method became known as the Average 
Density Method. 

Poettmann [8] derived an expression for calculating the sand- 
face pressure of flowing gas wells in which the variation with pressure 
of the compressibility factor of the gas is taken into consideration. 

Cullender and Smith [9] developed a workable procedure for 
calculating pressures in gas wells and pipelines which makes no 
assumptionsregarding either temperature or compressibility. Since 
no assumptions are made regarding the temperature and compressibility, 
this method is believed to give a more accurate prediction of the 
bottom-hole pressure. 

A method presented by Fowler [10] assumes a constant average 


temperature and involves integrated values of the gas law deviation 
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factor with pressure. It is a direct method of calculating static 
bottom-hole pressure. Sukkar and Cornell [11] extended Fowler's and 
Poettmann's methods and presented a general approach for calculating 
both static and flowing bottom-hole pressures for natural gas wells. 
However, Sukkar and Cornell only tabulated integral values for P, 
from 1.0 to 12.0 and T,, from 1.5 to 1.7 and thus the method has 
limited applicability. 

Dranchuk and McFarland [12] presented a technique for numerically 
integrating the differential form of the force-momentum balance in 
order to calculate pressure profiles in flowing gas wells for which 
the temperature profiles are known. The technique is similar to that 
presented by Young [13] except that it (a) employs a variable step- 
size, which reduces the number of steps required while increasing 
accuracy; (b) allows for the variation of both compressibility and 
temperature with depth in the kinetic energy term; (c) effects inte- 
gration in such a manner that both production and injection cases may 


be treated with one equation which has no singularities. 
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CHAPTER III 


THEORY 


All of the methods studied in this work are directly or 
indirectly derived from the general flow Equation (1) by the 
introduction of certain assumptions. The derivation of each method 
may be outlined as follows: starting with the general flow equation 


which may be written as: 
du 2 
Vdp + 5s + 2 dQ + + We = 1 
J. Vap f han Ly ay (1) 


One may replace V by 1/p, which reduces the equation to 


Since there is no shaft work, W.= 0; usually Au is small, 
consequently methods 1 to 5 assume that the kinetic energy term can 
be taken as zero. By assuming g = gc numerically, Equation (2) can 


be rewritten as: 


pes ig de + &, = 0 (3) 


i (© 


The lost work term, &,, can be estimated by the following 


equation: 


i 7 (4) 
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The modified gas law states that 


uy Ps ape W 
EN 8 @ZnRT was yw RT (5) 


Since density, po is defined as w/V", one can write 


Peet 
fi ObAZRT (6) 


This equation can be rearranged to yield 


The gas gravity is defined as 


eee vases cases (8) 
Mair». 28-96 


By substituting this equation and the value of R, the universal 


gas constant, into Equation (7), one obtains 
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Equations (4) and (9) can then be substituted into Equation (3) 


to yield 
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Before integrating, it should be noted that Q, and u are 
interrelated through the expression 
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where 


= 
oO 
ul 


flow rate in MMSCF/D 


‘= 
iT] 


velocity in ft/sec. 


Assuming Zo = 1 and substituting Equation (11) into Equation 


(10), one obtains 
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This equation is the starting point for the derivation of 
methods 1 to 5. The difference between these methods arises from 
the manner in which the temperature and the compressibility factor 
are handled so as to effect integration. These differences may be 


Summarized as follows. 


1. Average Temperature and Average Compressibility Factor Method 


Assumptions involved in this method are: 

(a) Steady-state flow 

(b) Single-phase gas flow 

(c) Change in kinetic energy is small and 
may be neglected 

(d) Constant temperature at some average value 

(e) Constant compressibility at some average 
value 

(f) Constant friction factor over the length 


of the conduit 


not 7 
- 
' 


to notiaviwst 4dy io 3nsoq corres ae 27 rt sae ‘ete 
mov? 2szrmw ebodisn e293 fisiatod sansisit : 
vote) yititoFzesyqmes wit baw cis crac 


ed Yer esonessttrh) s2ari MOTI HBSS He joa# 


burial so1964 yiiitdtezwrgn) sieiowh beg bing! 6408 some 
3315 borlirem 24 iid nf’ Bp we 
warts hte 

laden 

ons heli rs qor8n9.. ai tants ) 


Define 


3.04 
as BP (13) 


and 


f Qo <T> <Z> Po 


2 
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Substituting T = <T>, Z = <Z> and Equations (13) and (14) into 


(12), one obtains 


Multiplying Equation (15) throughout by a yields: 
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This is an ordinary differential equation and can be solved for 


the following boundary conditions: 


k = ee = Q 3 P Py 


k= Lo = Cal ; Pp = P. 


The solution may be written as: 
-2L -2L 
ee ple’ =i (17) 
val 
By multiplying Equation (17) throughout by e@ yields: 
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By defining 


ae aie 2GL 
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and substituting b and c into Equation (18), one obtains 


fQ2 <T>2 <JZ>2 p2 
e Ps : Ps Oe (elie e°) (19) 
29,D (0.0216) ites 
If one lets 
i = 14.65 psia 
Te =" 60°F = 9"9520°R 
D = d/123; where d = diameter in inches 


and substitutes into Equation (19), one obtains the working equation 
for the Average Temperature and Average Compressibility Factor Method, 
which is, 


0.5 
p Cae 25 QéG <T> <Z> fl (e© - 1) 
= e ee 
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The method consists of picking an initial guess for P,, 
calculating <Z>, at the average temperature and pressure, and the 
friction factor and calculating P, using equation (20). The result 
is compared with the initial guessed value. The procedure is 
repeated until a match in the value of P, is obtained. Data required 
for this method are gas gravity, flow rate, well depth, flow pipe 


diameter, wellhead pressure, wellhead temperature and bottom-hole 


temperature. 
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2. Average Density Method 


Assumptions involved in this method are: 

(a) Steady-state flow 

(b) Single-phase gas flow 

(c) Change in kinetic energy is small and 
may be neglected 

(d) Gas density is constant at an average value 

ROS On tea) i 2 

(e) Friction factor is constant over the length 

of the conduit 


(f) Velocity of fluid is constant at an average value 
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The average density can be written as: 
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By substituting Equation (23) into Equation (24), one 


obtains 
ye 4 x 10° QGP, 2 x 53.34 1 
(24 x 3600)nD? 93.34 Ty G Ps P. 
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Integration of this equation and substitution for T, yields 
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This is the working equation for the Average Density Method. 


The method of solution is similar to that described in Method 1. 
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3. Cullender and Smith Method 
Assumptions used in this method are: 
(a) Steady-state flow 
(b) Single-phase gas flow 
(c) Change in kinetic energy is small and 
may be neglected 
(d) Temperature gradient is known 
(e) Constant friction factor over the length 


of the conduit 


With the above assumptions, Equations (9) and (11) can be 


substituted into Equation (21) to yield 
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Multiplication of the right hand side of this equation by 


Pare ha 
(+5) / (75) produces 


Glo ia [rich hint tc OI stihes sto SRE: (30) 


2 P \2 H T | QoPo j 
TZ 2g-D° \0.0216 mT, 


bite iol 2 vevane os fan Lar. 
poh! 


Awint! ate hers — if (bh 7 


KJansT ond eve 705 nt sekcaks n or 7 
ne bo 


ad neo (11) ons (@) anofdeups sttcnicaaa ait 1s) 7 


1 
\ at 


bISiy. ot (TS) sorsanpa aint sstosrede 
ere *@ ay 


tS otn0 ST. BeLER 
- 4b (55) (| ee ot Sb + Se 
= 8 4) (arsine | ap ie Oe eae 


CS) 


: rt Pons ! 
(e cl eel i 
HY (natn a 0 ‘ 
| sys 
: _ Zi: : 
xd noted Hla a te 
Sy : a oth 
ie 
7 aT 7 


14 


Substituting 
Ge oe = Peale 
D = d/l2 
Po = 14.65 psia 
and i =e OcUEh 


dP (31) 


For convenience, one may define 


AWE 
Reale £93 
(5 + 666.6 Pi 


in which case Equation (31) can be written as 


sae Sale (32) 


This is the working equation for the Cullender-Smith Method. 
An evaluation of the integral over defined limits can be accomplished 
by numerical means. The method is rather tedious if a large number 
of increments is chosen for L. However, the authors suggest that by 
means of a two-step calculation and the application of Simpson's 
Rule, reasonable accuracy can be obtained. Thus equation (32) can 


be evaluated numerically using the following form: 


eGL_ eAP +4] +1) (33) 


an 
5 
7 ¢ 
ty 


Lin é : 
ane 


2 


(SE) 


botisalt dp tin-robest (02) ai jot vantage ah 
banat fqn0do, aa nee ental: aot ito 


vant 


varimusy spv8) ‘Bat btiepies ids neh — 3: 
a ‘get seancde 2% ona ue ot 


ie _ 
“et ; (2 40) ‘ oa ine tae 
a | cy ps ) i au T. 
7, ni b ges A : td 
pa Jee =) 
- - 


oe - = 
1D OnE ah, iG t-te t & 
6 yn id oi 
> : . 


15 


where 
P = Pes - Be 
and 
igs = wellhead pressure 
P, = bottom-hole pressure 


The method of solution consists of a trial and error evaluation 
of i I. and besa The bottom-hole pressure can then be calculated 


from Equation (33). 


4. Poettmann Method 
This method employs the following assumptions to derive 
practical equations from the basic energy balance: 
(a) Steady state flow 
(b) Single-phase gas flow 
(c) Change in kinetic energy is small and 
neglected 
(d) Temperature is constant at some average value 


(e) Friction factor is constant over the length 


of pipe 
(f) Velocity is constant at an integrated average value 
Ee a 
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Equation (21) can be written in terms of <u> 
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piles ‘ f<u>? 
Jdon | taal f+ Fes |e = 0 (35) 


Combining Equations (11) and (34) and substituting T = <T>, 


one obtains 
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in which case Equation (36) can be written as 
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Recall that 
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By substituting Equations (38) and (39) into equation (35), 


One can integrate the resulting equation to obtain 
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and substituting the values of T, and P, into Equation (40), one 
obtains the working equation for the Poettmann Method 
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Method of solution involves trial and error with procedures 


as described in Method 1. 
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5. Sukkar and Cornell Method 
This method uses the same assumptions as the Cullender and 
Smith Method except it makes the additional assumption that 
2 ==cons tant. 


Thus, by substituting T = <T>, Equation (31) can be written 


as 
GL 1 Z/P 
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d P 
Substituting Equation (37) into the above equation, one 
obtains 
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For convenience, one may define 
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in which case Equation (43) may be written as 
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If one defines 


P 
I = Vy Z/Py. dP 
SC1 02 Te r 
1+ BS 
r 
and 
I = Pre salt 
a f dP 
1 + al 
‘f 
where 
Dey = Sukkar-Cornell integral at point 1 
and 
we i Sukkar-Cornell integral at point 2 


then Equation (44) can be written as 
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This is the working equation for the Sukkar-Cornell Method. 
This is a direct calculation method which does not require trial and 


error. Method of solution consists of calculating ize from Equation (45) 
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and obtaining the value of P from tabulated data that will give 


this integral value. 
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6. Dranchuk and McFarland Method 
This method employs a technique to integrate the force-momentum 

balance equation numerically in order to calculate pressure 
profiles in flowing gas wells for which the temperature gradient 
is known. The method employs a variable step-size, which reduces 
the number of steps required while increasing accuracy. This 
method does not delete the udu term as do the other methods. 
It also allows for the variation of both compressibility and 
temperature with depth in the kinetic energy term. 


The working equation for this method is 
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By choosing the usual reference conditions, field units and 


applying the friction factor correlation 
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as suggested by Cullender and Smith [14], the appropriate values for 


the three coefficients become 
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and 
Cam = O.0L87489G 
where 
G = gas gravity 


Equation (46) along with the indicated values of the coefficients 
may be integrated to yield pressure profiles for both producing and 
injection wells. 

Generally, analytical solutions are preferred if possible to 
furnish a general solution to a problem. However, many phenomena, 
such as the one encountered in this work, are non-linear and the present 
mathematical methods proved to be well-nigh unworkable. Or perhaps 
analytical solutions are available but in such a form that is incon- 
venient for direct interpretation numerically. 

For the six methods studied in this work, closed-form solution 

available only for methods 1 and 2 but not for the other four 
methods. Therefore analytical solutions were developed for these two 


cases. For the Average Temperature and Average Compressibility Factor 
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Method, the starting point of the derivation is Equation (20). By 


differentiating this equation with respect to each variable studied, 


the following equations results. Adetailed derivation is shown in 


Appendix C. 


dP, 


2GL 
9 103.34 <T> ne en GR 12 Qo <T> <Z> ete 
ee i pee: qa 
(a Ni 
(a7) 
| 2GL | 
537 34e<T>: <7> 


25 GQ,fl <T> <Z> [e 


é 1] 


(48) 
Py 
: 2GL 
se Sy OLS te] fen of le (49) 
1 
kK, 
ee <> <I> -ki kp 2 “Ki ke [2@<Z> 
dT, QeTSeeTS* EIS PAR pe p 
r 
+ Ks os <T>2 2<Z> + Ks eps2eTs 
k 2417 | 5 ko 
2 P 
Yr 
a sh <I>*<> - — <Z> 
k, [9273 k3 uel 
nif alial zara |. ba eal eared le 
lf Y 
Ks (50) 


. | 
re sstbute sfdinbaew ml 29 ath ase ovo ps STP Ent in 


“ni mips et rotget96 seat ny ait a, i 


‘a 16 Ms 7 : ar 
; as AE) — nin st a : MO) A P pig icial a 


sh af 
ve 

7" et aoanA 

5 _ } "a J ; tl ad 


ir - 
(38) Whe. 
z ae 
a 
49) Ly 
= - 
m | sarah RTE ae 
1?, oF ‘ sie 7) 
: » A ' 
: ’ 
j 4] > 7 \ 7 : 
eO8 |S pte ehh Vass Ts | a 
A ea li ne ise is] | : | 


7 xs 


ep 


ES ei ata ib a ae 
an te a : >) 
an a ae 


7 


+a 


aa 
aif 


; Si 
es A 


23 


where 
i Fe; 
ye Aeal 
2 53734 
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Similarly for the Average Density Method, by differentiating 
P in Equation (27) with respect to each variable of interest, the 


following equations result. Adetailed derivation is shown in 


Appendix C. 
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Equations (47) to (57) can be used directly to calculate the 
effect of changing each variable in the calculated bottom hole pressure 
using either the Average Temperature and Average Compressibility Factor 


Method or the Average Density Method. 
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CHAPTER IV 


CALCULATION METHODS 


Computer programs were written for each of the six methods 
Studied in this work. Correlations for pseudocritical properties, 
compressibility factor, viscosity and friction factor of the flowing 
fluid are given in the subprograms. Appendix D gives the listings 
of all the computer programs and subprograms used in this study. 

Pseudocritical properties of the flowing fluid are calculated 
using either the fluid compositions or gas gravity. If fluid 
composition is available, Kay's mixing rule [15] is used to calculate 
the pseudocritical properties. The presence of sour gases are 
accounted for by using Wichert and Aziz [16] correction to the 
pseudocritical temperature and pressure. However, when the fluid 
composition is not available and where the gas gravity is the only 
known gas property, the pseudocriticals may be estimated from the 
correlation by Brown et al. [17]. Since this correlation is presented 
in graphical form, algebraic equations are written to fit the points 
on the graph to permit computer applications. If the mole fractions 
of the sour gases present are known, the pseudocriticals of sour 
natural gases may be estimated from the correlation by Carr et al. 
[18]. The estimated pseudocriticals are then corrected for the 
presence of C0, and H.S using Wichert and Aziz [19] method. 

The compressibility factor of dry sweet natural gas is 


calculated by using the reduced form of the B-W-R equation of state 
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as suggested by Dranchuk, Purvis and Robinson [20]. This method can 
be used for sour natural gases by using modified pseudocritical tem- 
perature and pressure, as suggested by Wichert and Aziz [21]. This 
correlation was tested and extended to predict compressibility factors 
of natural gases up to a reduced pressure of 30. 

Fluid viscosity is calculated using the method suggested by 
Abou-Kassem [22]. 

Friction factor is calculated by the method suggested by 
Dranchuk and Abou-Kassem [23]. 

A numerical scheme was set up to evaluate both the Poettmann 
and Sukkar-Cornell integrals. The integration method employed 
Simpson's Rule and Newton's 3/8 Rule as is decribed in Appendix B. 

For the Average Temperature and Average Compressibility Factor 
Method, Equations (47) to (57) can be used to study the sensitivity 
of each parameter. For example, by substituting the values of G, 


<T>, <Z>, f and L into Equation (48), the following equation results: 


1 Qo 
—— = 6990.346 —— (58) 
dQo Py 


However, the assumption has to be made that f is constant within the 
range of Qy studied. Equation (58) can be integrated to give a general 
solution if the boundary conditions are known. However, for sensitivity 
analysis purposes, Equation (58) can be integrated within the range of 


limits of interest, say 10% to give, 


Pe - PL = 6990.346(Q5 - Q) (59) 


Vd Haszobou2 Sinitom ant yd Bagetuals9 re 


anamsteod ect dio? sieulsve of . 7: 2 ait iagite 


bayotqms boryan scab i ne Be 


8 xthneagA of badinrogb: zt be Stun B\E 2 ‘Hoswail ne afui peak 
want “rT tdkezavqmod adedsvA tan hs A ct \ 38 gers 


WiVistamse add \uite vi be2u 9d Tso (yay oy pe j | 
,o Ta 2aUTsy ‘sas OM Sad Pl 2uue va so Tqnexs ah on re i 


:2Sfveay notssups patwoligi Baz , (ah) pil phos? esi> ** 


(82) 


=) 


a 


28 


where 


Equation (59) can be rearranged to give, 


6990.346(Q, + Q,)(Q, - Q,) 


AP : (60) 
P, + P, 


where 


Mas wate P; 


The above procedure can be applied to Equations (47) to (57) 
to study the sensitivity of each parameter on the calculated bottom- 
hole pressure. However, the assumption that changing the parameter 
of interest only affects the calculated P,;, has to be made in order 
that integration can be performed. For equations involving the gradi- 


ent af and a, graphical methods were used. 
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CHAPTER V 


RESULTS AND DISCUSSIONS 


Two flowing gas wells with measured bottom-hole pressures were 
used in this work for the purpose of study and comparison. Each 
variable was varied incrementally up to +10% and the bottom-hole 
pressure was calculated using all the six methods. The well data 
are shown in Table 1. These well data and an absolute roughness of 
0.0006 inch were used as the basis to calculate the deviation in the 
calculated bottom-hole pressure upon varying the input variable. 

Tables 2 and 3 show the comparison of the bottom-hole pressure 
calculation methods for Wells 1 and 2 respectively. For Well 1, all 
the methods give bottom-hole pressures less than the measured value. 
Differences are within the range of 32 to 42 psia. For Well 2, all 
the methods give higher pressures than the measured value except for 
the Poettmann method which gives a bottom-hole pressure of 55 psia 
less than the measured value. Differences are within the range of 
20-30 psia for the other methods. 

The sensitivity analysis of the various methods to variations 
in input variables are given in Tables 4 to 15 and the results are 
plotted on Figures 1 to 20 (Appendix A). Table 16 gives the summary 
of the maximum percent deviation in the calculated bottom-hole 
pressure for the six methods for both wells. Table I gives the degree 
of sensitivity (in increasing order) and the associated maximum percent 


deviation in the calculated bottom-hole pressure of the various methods 
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TABLE 1 


SUMMARY OF WELL DATA 


Well depth, ft. 


Flow rate, MSCF/D 


Tubing 19D... sin. 


Gas gravity 


Wellhead temperature, °F 


Bottom-hole temperature, °F 


Temperature gradient, °F/ft. 


Measured B.H.P., psia 


Mole % CO2 


Mole % No 


Mole % H2S 


68.0 


115.0 


0.0099 


0.00 


WELL 2 


11029.0 


ASSEN IO 


22992 


0.702 


170.0 


243.0 


0.0066 


3690.4 


2.84 
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DEGREE OF SENSITIVITY AND THE ASSOCIATED MAXIMUM PERCENT 


DEVIATION IN THE CALCULATED BOTTOM-HOLE PRESSURE OF THE VARIOUS 
METHODS TO DIFFERENT INPUT VARIABLES FOR WELLS 1 AND 2 


Dn oO FP W 


VARIATION IN DEPTH 

WELL 1 
Poettmann (1.63%) 
Sukkar and Cornell (1.69%) 
Dranchuk and McFarland (1.70%) 
Average Density (1.70%) 
Cullender and Smith (1.70%) 
Average T and Z (1.71%) 


VARIATION IN FLOW RATE 


WELL_1 
Dranchuk and McFarland (1.13%) 
Cullender and Smith (1.15%) 
Average Density (1.15%) 
Average T and Z (1.16%) 

Sukkar and Cornell (1.32%) 


Poettmann (1.36%) 


WELL 2 
Poettmann (3.04%) 

Dranchuk and McFarland (3.16%) 
Sukkar and Cornell (3.20%) 
Average Density (3.20%) 
Cullender and Smith (3.23%) 
Average T and Z (3.26%) 


WELL 2 
Dranchuk and McFarland (1.88%) 
Average Density (1.95%) 
Cullender and Smith (1.96%) 
Average T and Z (1.99%) 

Sukkar and Cornell (2.23%) 


Poettmann (2.51%) 
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TABLE I (CONTINUED) 


VARIATION IN WELLHEAD PRESSURE 


WELL_1 
Poettmann (8.57%) 

Sukkar and Cornell (8.82%) 
Average T and Z (9.02%) 
Cullender and Smith (9.03%) 
Average Density (9.03%) 
Dranchuk and McFarland (9.08%) 


WELL 2 
Poettmann (7.44%) 

Sukkar and Cornell (7.78%) 
Average T and Z (8.0%) 
Dranchuk and McFarland (8.03%) 
Cullender and Smith (8.04%) 
Average Density (8.05%) 


VARIATION IN BOTTOM-HOLE TEMPERATURE 


WELL 1 
Sukkar and Cornell (0.051%) 
Average Density (0.097%) 
Cullender and Smith (0.099%) 
Average T and Z (0.10%) 
Poettmann (0.145%) 
Dranchuk and McFarland (0.68%) 


WELL 2 
Sukkar and Cornell (0.31%) 
Average Density (0.39%) 
Cullender and Smith (0.40%) 
Average T and Z (0.41%) 
Poettmann (0.62%) 

Dranchuk and McFarland (1.38%) 
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TABLE I (CONTINUED) 


VARIATION IN WELLHEAD TEMPERATURE 


WELL_1 
Poettmann (0.056%) 

Average T and Z (0.058%) 
Cullender and Smith (0.059%) 
Average Density (0.061%) 
Sukkar and Cornell (0.063%) 
Dranchuk and McFarland (0.65%) 


VARIATION IN GAS GRAVITY 
WELL_1 
Poettmann (1.58%) 
Sukkar and Cornell (1.91%) 
Average Density (1.95%) 
Average T and Z (1.97%) 
Dranchuk and McFarland (1.99%) 
Cullender and Smith (2.34%) 


WELL 2 
Sukkar and Cornell (0.24%) 
Poettmann (0.25%) 

Average T and Z (0.28%) 
Cullender and Smith (0.29%) 
Average Density (0.31%) 
Dranchuk and McFarland (1.36%) 


WELL 2 
Poettmann (2.88%) 

Dranchuk and McFarland (3.50%) 
Sukkar and Cornell (3.71%) 
Average Density (3.93%) 
Cullender and Smith (3.96%) 
Average T and Z (4.0%) 
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ORDER OF SENSITIVITY OF EACH VARIABLE AND THE ASSOCIATED 


MAXIMUM PERCENT DEVIATION IN THE CALCULATED BOTTOM-HOLE PRESSURE 
FOR THE SIX METHODS 


(1) 


AVERAGE T AND Z METHOD 

WELL 1 
Wellhead Temperature (0.058%) 
Bottom-hole Temperature (0.1%) 
Flow Rate (1.16%) 
Depth (1.71%) 
Gas Gravity (1.97%) 


Wellhead Pressure (9.02%) 


AVERAGE DENSITY METHOD 


WELL 1 
Wellhead Temperature (0.061%) 
Bottom-hole Temperature (0.097%) 
Flow Rate (1.15%) 
Depth (1.70%) 
Gas Gravity (1.95%) 
Wellhead Pressure (9.03%) 


WELL 2 
Wellhead Temperature (0.28%) 
Bottom-hole Temperature (0.41%) 
Flow Rate (1.99%) 

Depth (3.26%) 

Gas Gravity (4.0%) 


Wellhead Pressure (8.0%) 


WELL 2 
Wellhead Temperature (0.31%) 
Bottom-hole Temperature (0.39%) 
Flow Rate (1.95%) 

Depth (3.20%) 
Gas Gravity (3.59%) 


Wellhead Pressure (8.05%) 
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TABLE II (CONTINUED) 


CULLENDER AND SMITH METHOD 
WELL 1 

Wellhead Temperature (0.06%) 

Bottom-hole Temperature (0.10%) 

Flow Rate (1.15%) 

Depth (1.70%) 

Gas Gravity (2.34%) 

Wellhead Pressure (9.03%) 


POETTMANN METHOD 


WELL 1 
Wellhead Temperature (0.06%) 


Bottom-hole Temperature (0.145%) 


Flow Rate (1.44%) 

Gas Gravity (1.58%) 

Depth (1.63%) 

Wellhead Pressure (8.57%) 


WELL 2 
Wellhead Temperature (0.29%) 
Bottom-hole Temperature (0.40%) 
Flow Rate (1.96%) 

Depth (3.23%) 
Gas Gravity (3.96%) 
Wellhead Pressure (8.04%) 


WELL 2 
Wellhead Temperature (0.25%) 
Bottom-hole Temperature (0.62%) 
Flow Rate (2.51%) 

Gas Gravity (2.88%) 

Depth (3.04%) 

Wellhead PHECeLEe (7.44%) 
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TABLE II (CONTINUED) 


SUKKAR AND CORNELL METHOD 


WELL 1 
Bottom-hole Temperature (0.05%) 
Wellhead Temperature (0.06%) 
Flow Rate (1.32%) 

Depth (1.69%) 

Gas Gravity (1.91%) 

Wellhead Pressure (8.82%) 


DRANCHUK AND McFARLAND METHOD 


WELL_1 
Wellhead Temperature (0.65%) 
Bottom-hole Temperature (0.68%) 
Flow Rate (1.13%) 

Depth (1.70%) 

Gas Gravity (1.99%) 

Wellhead Pressure (9.08%) 


WELL 2 
Wellhead Temperature (0.24%) 
Bottom-hole Temperature (0.31%) 
Flow Rate (2.23%) 

Depth (3.20%) 

Gas Gravity (3.71%) 

Wellhead Pressure (7.78%) 


WELL 2 
Wellhead Temperature (1.36%) 
Bottom-hole Temperature (1.38%) 
Flow Rate (1.88%) 

Depth (3.50%) 

Gas Gravity (3.50%) 

Wellhead Pressure (8.03%) 
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to different input variables for Wells 1 and 2. Table II gives the 
order of sensitivity (in increasing order) of each variable and the 
associated maximum percent deviation in the calculated bottom-hole 
pressure for the six methods. 

It can be seen that for both wells, the calculated bottom-hole 
pressure is most sensitive to variation in wellhead pressure and 
least sensitive to temperature variations. However, Dranchuk and 
McFarland Method is very sensitive to temperature variations (about 
10 fold for Well 1 and 5 fold for Well 2) when compared to the other 
methods. Among the six methods studied, Poettmann Method is least 
sensitive to variations in depth, wellhead pressure, wellhead temper- 
ature and gas gravity for Well 1. The same is true for Well 2 except 
for the variation in wellhead temperature. Methods 1, 2, 3 and 6 
show the same trend for both wells regarding the order of sensitivity 
of the variables. They are, in the order of increasing sensitivity, 
wellhead temperature, bottom-hole temperature, flow rate, depth, gas 
gravity and wellhead pressure. 

Poettmann Method and Sukkar-Cornell Method also follow the 
same order for Well 2. For Well 1, gas gravity is more sensitive 
than depth in the Poettmann Method and bottom-hole temperature is 
more sensitive than wellhead temperature in the Sukkar-Cornell Method. 
In all the six methods, the sensitivity of each variable increases, as 
the depth and flow rate increases, except for wellhead pressure. 

Tables 17 and 18 show the effect of absolute roughness on the 
calculated bottom-hole pressure for the four methods studied and the 


results were plotted on Figures 19 and 20. The results show that the 
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Poettmann Method is more sensitive to variation in absolute roughness 
than the other three methods studied. Thus if the absolute roughness 
is in error, an error will be introduced in the calculated friction 
factor and thus the calculated bottom-hole pressure. This may explain 
the inaccuracy of this method when compared with the other methods. 

In order to study the effect of using linear interpolation of 
tabulated integral values for the Poettmann Method and the Sukkar- 
Cornell Method, the bottom-hole pressure was hand calculated for 
Well 1 using the two methods. The integral values for both methods 
were estimated by linear interpolation from the table of integral values 
established in this work. The comparison was made for +1.0% and +102% 
variation of each variable and the results are shown in Tables 19 
and 20. 

Table 19 shows that the error introduced by using linear inter- 
polation in Poettmann Method is small. The base case gives an error 
of only 0.006%. This is to be expected since the step-size taken in 
P, and Ty is small (0.05) in evaluating the integral. However, this 
is not the case in the Sukkar-Cornell Method. The error introduced 
by using linear interpolation in the base case is 6.55% which is very 
significant. The reason for this is that the step-size taken in B is 
large (5.0). Also in the Sukkar-Cornell Method, the bottom-hole 
pressure is calculated from the product Py times Pc. where P, is 
obtained by linear interpolation. Since P, is of the order of 600 to 
J00=DS ta, even a small error in the calculated P, will introduce 
a significant error in the calculated bottom-hole pressure. Reducing 


the step-size in B will require too much computer time and becomes 


uneconomical. 
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The effect of 10% change in input variable on the calculated 
bottom-hole pressure using the equations derived for Methods 1 and 2 
are summarized in Table 21. It can be seen that the numerical results 
are in good agreement with that obtained by analytical methods. 

The sensitivity analysis range of +10% studied in this work is 
sufficiently large and should cover the likely measurement error that 
would normally occur in practice. The well depth is measured very 
accurately during well completions. The error introduced here should 
be minimum and should be in the order of less than 0.05% in a 10,000 
ft. well. The tubing head pressure is normally measured by a dead 
weight gauge or Barton gauge. The accuracy of these devices should be 
within 2% if properly maintained and calibrated. Temperature measure- 
ment is normally measured at surface with thermocouples or temperature 
gauge and with bottom-hole temperature bombs at the sand face. These 
measuring devices are usually accurate if properly calibrated and the 
error introduced should be less than 2%. Flow rates are commonly 
measured with an orifice meter. If properly designed and located in 
the proper location to avoid turbulence, accuracy should be within 
1 to 2%. Gas gravity is usually calculated from the gas composition 
obtained from 9JaSachromatograph. If the chromatograph is well main- 
tained, the analysis should be very accurate and the calculated gas 
gravity should be within 1%. 

In summary, the instrumental error introduced in measuring the 
various parameters is generally less than 2% if the equipment jis 
properly maintained and calibrated. The biggest error that normally 


occurs in field operation is in fact human error. Therefore, if the 
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TABLE 21 


COMPARISON OF RESULTS 


NUMERICAL SOLUTION VS ANALYTICAL SOLUTION 
USING DATA FROM WELL 1 


<T> AND <Z> AVERAGE DENSITY 
METHOD METHOD 
VARIABLE NUMERICAL ANALYTICAL NUMERICAL ANALYTICAL 
(10% CHANGE ) RESULTS RESULTS RESULTS RESULTS 
DEPTH (L) ZA 1.70 sea) 1.68 
FLOW RATE (Q,) 1316 WZ TNS 1.14 
WELLHEAD P. (P2) 9.02 8.87 9203 8.61 
B.H. TEMP. (11) -0.092 -0.081 -0.086 -0.083 
W.H. TEMP. (T2) -0.055 -0.048 -0.057 -0.053 


GRAVITY (G) leo7 1.70 os ey 
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calculated bottom-hole pressure calculated by any one of the six 
methods is significantly different from that of the measured value, 
the field engineer should double check the various input parameters, 
especially the wellhead pressure and flow rate, which are most sus- 


ceptible to human error, to ensure that the data is correct. 


s tudy: 


CHAPTER VI 


CONCLUSIONS 


The following conclusions are made as a result of this 


Calculated bottom-hole pressures by all the methods 
Studied are most sensitive to variations in wellhead 
pressure and least sensitive to variations in 
temperature. 

The percent change in calculated bottom-hole pressure 
shows the same trend for all the methods except for 
variations in bottom-hole and wellhead temperature 

in the Poettmann Method. 

The Poettmann Method is more sensitive to variations 
in absolute roughness than the other three methods 
studied. 

At high flow rates, the kinetic energy term may be 
significant, thus the Dranchuk and McFarland Method 
should give the best results. 

Linear interpolation can be used with confidence in 
the Poettmann Method. 

Linear interpolation should not be used in the Sukkar 


and Cornell Method. 
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For the Average Temperature and Average Compressi- 
bility Factor Method and the Average Density Method, 
the numerical results obtained are in good agree- 
ment with the analytical results using the analyt- 


ical solutions developed in this work. 
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CHAPTER VII 


RECOMMENDATIONS 


The following recommendations are made: 


Collect flow data for more wells with different depths 
and flow rates and perform bottom-hole pressure calcu- 
lations with the six methods studied. This will enable 
one to make a better comparison of the different methods 


and the limits of their applicabilities. 


Investigate possible instrumental error involved in 
various measuring devices that are used to measure the 
input variables required to do bottom-hole pressure 


calculations. 
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NOMENCLATURE 


inside pipe diameter, inches. 
inside pipe diameter, feet. 
Moody friction factor = 4 f', dimensionless. 


Fanning friction factor, dimensionless. 


local gravitational acceleration, 32.174 ft/sec’. 


conversion factor, 32.174 1b/ft/1besec*. 

gas gravity (air = 1). 

length variable, feet. 

lost work, ft.1b-/Ib... 

length of flow string, feet. 

molecular weight, 1b /1b-mole. 

pressure variable, psia. 

pseudocritical pressure, psia. 

pressure at reference conditions, 14.65 psia. 
pseudoreduced pressure, dimensionless. 

flow rate, MMSCF/Day. 

Universal gas constant, 1545 ft. 1b¢/Ib-mole°R 
temperature, °R 

pseudocritical temperature, °R 

temperature at reference conditions, 520°R 
average temperature, °R 

velocity, ft/sec. 

average velocity, ft/sec. 


specific volume, ft°*/1b,. 
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volume, ft?. 
mass velocity, Ib /ft*-sec. 
mass, 1b, 
Shaft work, ft-1b,/1b : 
m 
elevation 
compressibility factor, dimensionless. 
compressibility factor at standard conditions, 
dimensionless. 


average compressibility factor, dimensionless. 


gas density, Ib /ft*. 
gas density at standard conditions, Ib /ft®. 


average gas density, 1b /ft®. 
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APPENDIX B 


This integration scheme is listed in the computer program as 
Subroutine QSF in Appendix D. This subroutine performs the integration 
of an equidistantly tabulated function by Simpson's rule. To compute 


the vector of integral values 


Xj 
Le = Z(x;) = ff y(x) dx 
a 
with iges [R22 saan) 
Sd (i-1)h 
for a table of function values y; (ies, yee ee Ny ediVven atvequi= 
distant points 4 ibaa Ciel )ns (1) =" 15 2a een) FP eS mpSOnas: rule 


together with Newton's 3/8 rule or a combination of these rules is 
used. Local truncation error is of the order of h® in all cases with 
more than three points in the given table. Only Z, has a truncation 
error of the order h* if there are only three points in the given 
table. No action takes place if the table consists of less than three 
sample points. 

The function is assumed continuous and differentiable (three 
or four times, depending on the rule used). 

Formulas used in this subroutine (Z; are integral values, y; 


are function values) are: 
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: as fan (1.25 Vay + 2y; 0.25 Ys) Reet hs re, ah) 
ef h . 

z. = Zi» “2 3 (V5.2 + ay 54 + y;) Simpson's Rule. . . (2) 
: 3 

. Zi 3 + oh (¥5_3 + 3Y 5.2 1 3Y 5) + ¥;) fa ee (3) 


Newton's 3/8 Rule 


h 
5 A ea : : ; ? } 
Z Z 3 (Vj_5 + 3-875 y;_4 + 2.625 y. 4 + 2.625 y,_, 


Combination of (2) and (3) 


Sometimes formula (2) is used in the following form 


Local truncation errors of formulas (1) - (4) are respectively: 


] UU) a 
Ree = 24 h*y (Eq) [Ei in(x, 7. X41) 
] soad 7 
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APPENDIX C 


For the Average Temperature and Average Compressibility Factor 
Method and the Average Density Method, closed-form solution is avail- 
able. Therefore analytical solutions can be obtained for these two 
cases. Generally, analytical solution is preferred if possible unless 
the final equation becomes too complicated and unmanageable. 

This appendix gives the detailed derivations of the analytical 
solutions for the two cases above. The final equations can then be 


used to check the accuracy of the results obtained by numerical methods. 


Average Temperature and Average Compressibility Factor Method 


1. Effect of well depth, L. 


The working equation is: 


24.99 Q2 G<T><Z>fL(e° - 1) 1? 
p) = e& pe +} 98 (1) 


dec 


: 2GL 
Se RS eee 


Let k, ps 


‘ 26 
KD = B98 gdel><2> 
25 Q)<T><Z>GF 


5 


d 


and eS 


Equation (1) can be written as follows, 


zzatny efthazog ty ligase @ 
| Sen © dade is (sh 

VG lone ony To gouttevted Nel isto ade zovip he 
ov isc) feo enorfeyns Fant? edT  /ovnds teeam gabe 
-aborttan. [astremn ve tantnado etfnaey oad a _ 


7 7 
a) 


bons Auitasd_yitnh IWtegenginay, 
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By differentiating Equation (2) with respect to L, one obtains, 


koL KoL 
dP, = kK, k, e : fi k, e : (3) 
dL 2P 


By substituting the values of k,, k,, and k, into Equation (3) 


one obtains 


2GL : : 
Le Mme "aati ESL APRY ES ee CR by 
dL oP, 
2. Effect of flow rate, Qo. 

Let k, Es ef Be 

d ae 25 GLf<T><Z>(e° - 1 
ty at disc 
where 


7 26L 
- 53. 34<T><Z> 


Substituting k, and k, into Equation (1) obtains 
L 
i eA cr ht k,Qo)* (5) 
By differentiating Equation (5) with respect to QW One obtains, 


at i SaMo i . (6) 
dQ (k, + k,Q°)? 


By substituting Equation (5) into Equation (6), one obtains, 
2GL 
53. 34<T><Z> 
Te<7> (e - | 
dP kQ, 25 GLFQ,<T><Z> ( ) 


aye = 2 _____§_— (7) 
dQ, P, Bs 


= aes r) 
rg ims 2 uh 
is 
¢ ms, 2 bee y 7 + 


a ey vol pe 
(8) nok ea cnt 8 808 eG a 


dip cai... eae he i 


ie 


hens 
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3. Effect of wellhead pressure, Po. 


25 Q2 G<T><Z>fL(e° - 1) 
fet kis —OoLe) 9 }-5 212 


2GL 
53. G4<le<L> 


Substitute k into Equation (1) to obtain 


pie (6° Po | (8) 


1 
By differentiating Equation (8) with respect to P,, one obtains, 


dos She nie 
—t = Wt. (9) 
dP» (e- PSF k,)? 


By substituting Equation (8) and value of c into Equation (9), 
One obtains 


2GL ) 


GE Apter cne borBd<iagze (10) 


dP, Py 


4. Effect of bottom-hole temperature, T, and wellhead temperature, To. 


fete Kie= aR 
oe 2GL 
oe Tee yesy! 
25 Oh GfL 
and k, lian a oa ee 
d 


Substituting k,, k2, and k; into Equation (1), one obtains, 
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(11) 


Equation (11) can be rearranged to yield 


Ke k, 
ak Pa pec <[e-</>" ao lathe - Ka eps2ez>2 (12) 
2 2 
Since 
Gis, = eae: 
2 
and 


dT, vA dT. 

therefore 

pole OSL de oie 
dba #7 &. d<T> edie 2 31> 


By differentiating Equation (12) with respect to T,, One obtains 


k, k, 
ee fale ee “Ki ky <ToeZs rk, ky see 
aT, ; aT, Sea d<Z>*<T> d<T> Pr. 
= Se 
+ Ka Sa aeee anlast> [ae + Pee as ><> 
Ko o<T> P,. 5 
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antedda eno .,1 03 tasqzar sew (81). OTS SEPS "QUT Cae 
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Equation (13) can be rearranged to yield 


kK, 
aP a gP - Benes 3K, kK, Keiko 1o5L2 
dT) dT. Q<Z><T>* — 9<Z><T>? | a<T>]} p 


K, 
_ + — <Z>*%<T> 
A Ke 


k 
eee <Z><T>? ee 
ko eka [es 


epee ea ae ee -} _ 
2 2 2 a p 


k 
= mee laa =] 


- Pa ee f(2e.) 
d<T> |} p k : 


Re 2 


Seg ei hect Of gas Gravity.ca. 


et uk = Pp: 
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=) gana le ty 
2 D334 L> 


25 Q6<T><Z>fL 


5 


d 


Substituting k,, kz, and k3 into Equation (1), one obtains, 


k k 
p iat k, ok2G 3 koG 8 (15) 


1 


By differentiating Equation (15) with respect to G, one obtains, 


dP, _ ef2M(ky ky + ks) | (16) 
dG Ae 


Substitute value of k,, kz, and k; to obtain 
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2GL ale E 25 Q5<T><Z>fL 
dP, ; ev3- 34<T><Z> ACCES pe Ep ores A 
dG Bis 
Average Density Method 
1. Effect of well depths (= 
The working equation is: 
P, + 9.374 x 107° GL a a 
Pp = + 9, x i ——- + 
: 3 3 T1Z) To2Zo 
Z ] 
, 24.99 Qo GLf (18) 
ds§ Py k PS 
Ty Ayes lad, 
Bene ee ey <0? le 
24.99 06 Gf 
aS eer cee 
Py P, 
and k3; = + 
: ; Us TZ, 
Substituting ki, k2, and k3 into Equation (18) yields 
Sap) Fete (19) 
Py a P, Ja 5>3 k 


3 


Equation (19) can be differentiated with respect to L to give 


dP, 
ee ie tore Un came ce ena (20) 
T, rye tei k2°T 575 dL 
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By grouping terms and substituting the value of ki, ko, and k3, 


Equation (20) can be rewritten as 


24.99 Q2 Gf 
05374 6. 10e 1G aan wie. Q 
TZ, T3525 5 P, P, 
dp d i 
NL PS Pe es 2 Set oe be OUI OL (21) 
dL 1 24.99 Q° GFL 
Pe | =}. 973744108) GE 
12,145 | Oa | 
es ToZ> 
PEt CeCtsOtme1OW fateemdnn 
Let k, = 9.374 x 107? GL 
, = 24.99 GLf 
es d° 
P P 
and ok, = —=+—+ 
TrZy- @laZe 
Substituting ki, k2, and k3 into Equation (18) yields 
k 2 
Pp, = P, +k, kz +— Xo (22) 


3 


Equation (22) can be differentiated with respect to %w to give 


dP dP 
k LUE ie 
dP, * dQ, ie k, k, 95 dQ 


(23) 
dQ, ey Koa k3 Z,T, 


By grouping terms and substituting the value of k,, kz and k; 


into Equation (23), one obtains 
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33 


Ot 


= 
49.98 9537) 4axOnewGe 
Py + P, 
ae = atten (Oeis Sea whiOe 1727 as eae ae (24) 
dQ, _ 9.374 x 1079 GL, 24.99 GLFQ3 
La P, 2) 2 
detZin ee 
Tada ToZo 


Effect of wellhead pressure, P,- 


let k, = 9.374 107? GL 
24.99 03 GLf 
4 = = 
d 
P P 
and k, = —— 
Ty Zy T2Z> 


Substituting k,, k,, and k,; into Equation (18) yields 


k 
Pames WP. Foky byt a (25) 
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Equation (25) can be differentiated with respect to P, to give 


dk k, dk 
Sines. Pray tae eas (26) 
dP dP» k; dP2 
Since 

a, 

dk dP ] P OZ 
_3 ‘= ea — [=| (27) 
dP, Pa Gp declare le 
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Equation (26) can be rewritten as 


k, \| dP, 1 Pana BZ 
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By grouping terms and substituting the value of k,, k,, and k3 


into Equation (28), one obtains 


| 24.99 2 GLE 
1 = = [9-374,3 TOP GL - valle (2 =, | 
Fi} P, 
anes rh Va 

shee 24.99 ae GLf 

1} - W027 10 ee 

12, n Po | 
T,4, 1,2, 


4. Effect of bottom-hole temperature, T; and wellhead temperature, T2. 


bet. ky = 9.3746 1092 Ge 
24.99 Q6 GLf 
k> Woe 
d° 
p P 
Tidy pe To Zo 


By differentiating k,; with respect to T,, one obtains 


dP, 
dk aT P 
nat = peal -=>> - (30) 
dT, ieee ees TZ 72 = 
1 
Substituting k,, kz, and k; into Equation (18) yields 
k, 


3 


Equation (31) can be differentiated with respect to T, to give 


ana + - = (32) 
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By substituting Equation (30) into Equation (32) and rearranging, 


one obtains 


al ba 


dP TZ Tae 7 eal 


dT 
ee sic 
TZ, K, 


Substituting the value of k,, k,, and k,; into Equation (33) 


yields 
24.99 Q* GLf ) ] Pee 
=} 9:374*« 10 * 7G = ae ene — + —|— 
d5 rae Pa TZ, qT, Zy oT, p 
dP M12) ToZe 1 
dT, : 24.99 yr GLf 
ie Set Gl ee 
T,2, d> | Py + P, | 
TZ; Tis Z5 
Pua Wee, Be (34) 
Similarly, 
2 
: ee S10 Ge | | & et | 
P P Tada lel Zot\8T 
dé Uae 2 DED 2 2 2! p 
dP _ Bent Lo < 
dT, . 24.99 Qe GLf 
1 (| 9.374 * 108 GL - are a aan) 
1,2, ToZ2 
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Dee ineGLEOt. das gravity. G:- 


Boteek, =§50-3740 10s aL 

24.99 Q2 Lf 

k, = ——_*— 

2 
d° 
p 
3 = 1 + eran 
TZ lees 


Substituting ki, k2, and ks into Equation (18) yields 


k 
P, = P, +k, ks Gt Ss (36) 


3 
Equation (36) can be differentiated with respect to G@ to give 
dP k,G dP, k, GaP, 


ee ak heat Sy ER (37) 
dG T, APACER KE T17, kk, dG 


wr 


Substituting the value of k,, k,, and k, into Equation (37) 


and rearranging yields 
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GY? 


MAIN PROGRAM 


J 


ate tes! U He sts sy! 4 J, 
PAE AE AE HEAT HE IS IIS IC HE ACC HE CI He He IC ke aK kk oie oc ole ake oe ake ste ale ste ole ste ole ae ote ate ok Safa endian bs 


COMPRESSIBILITY FACTOR METHOD, 
THIS PROGRAM CALCULATES THE EFFECT OF 
VARIATIONS IN THF INPUT PARAMETERS ON 
THF CALCULATFD BOTTOAM-HNLFE PRESSURE. 
X(K)=1 VARIATION IN PFPTH 
xk X(K)=2 VARIATION IN FLOW RATE 


oh ee ae 
vs a" WW 
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' 
ep erage a Pay reed 


sc 
* SFNSTITIVITY ANALYSIS OF AVERAGE TEMPERATURE AND *e 


s X(K)=3 VARIATION IN WFLLHEAD TURING PRESSURE * 
as X(K)=4 VARTATION IN BOTTOM-HOLE TEMPERATURE * 
*6 X(K)=5 VARIATION IN WELLHEAD TEMPERATURE *e 


A J A A, =! 
HEE HE ASHE OE IE DIS TC ETE DIC BKC OE IC aHe ole ate aie ole ae oie ole sis ois ate ak aie ale ste ak ale ole ste ote ate ate akc ale ate ale ake ak ale ole oe ote aioe ok ote of 


NOMENCLATURE 


YMOEF*=- MALE FRACTION 

YM = MOLECULAR WEIGHT 
TEME@R=PGCRITICAL TEMPERATURE. DEGREE R 
PRECT) H=ACRITICAL PRESSURE, PSTA 

OTAMM ASP INSIDE SDIAMFTER® DF MARLOW! PAPESINERES 
TSD = WFLLHEAD TEMPERATURE, DEGREF F 
TOD =TBOTTOM=HOREPTEMPER ATURE ,NEGREEFE 
PSN ="WELEHEADS PRESSURIES PSITA 

@) =" FUN RATE 

TRC =-RFNUCED TEMPFRATURE 

PRC REDO EPRRRE S SURE 


/TEMER= PSEWVDOCRITILCAL TWEMPERATURE. DEGREE R 


APREC =ePSEUDNCRITICAL PRESSURE, PSTA 
SG =F CASH GRAVITY 

NENSS se =O DENSTTY 

REN =F REYNOL DSONUMBIEE 

AVEMW = AVFERAGF MOLECULAR WEIGHT 

ViISMe FsFVISCOSTIN GP 

TEMC = CRITICAL TEMPFRATURF,NFGREE R 
PREECE. S=/CRITICARE PRESSURES PSI A 


NIMENSTON YMOLF(16),YMW(16),TFMC( 16) ,PRFEC(16),A(8) 


Ho ECON) 


MART AW A/0- 6150623 751.9467099 506 91832772950. 53530711, 
1 OG? 32032 , 001 04868144 0568 15:7:00150.68446549/ 


FED. OOO6 


REAM AND WRITE INPUT DATA 


WRITE (6,208) 
20 BP EORMAT{ ' 1") 


REAND( 5,201) (YMW(T),TFMC(1) sPRFEC(IT),1=1,16) 


201 FORMAT(3F10.5) 
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MAIN PROGRAM ese (COUNTED) 


WRITE (Gs 212) 


Ne. POURMATCL7 £119 29x, MOL suites te OX si CR 1 oe EME. ly 
Mer Int le PRES S59 
TATAS DEG. RANKING, 1 2X <u ST At /) 
WRITEUCEGs 202) ¢(YMW(T ) 5 TEMC(1) sPREC(I )»y1=1,16) 


292 FORMAT(10X,*CN2 OR 206599 LOXs IND 3 h20 eb 7 
10%, §H2S 
DpeP 20.59 7, 1 OX! Ol te 3E 20 65 of LOX Ce ',3F20.5, 


fe LOK, "ES 
1 Ms SE ZO D's / p LOXG* 1 =C4 13 3F 206.5975 10%s 'N=CS ' 
Ke SECU De op LOK, | f= 
1C5 ",3F20055/351 0X, 'N-C5 Wet 20 sD 97 aX pl Co ' 
*%y3F20.59/510X, 
1eN= Co Or ZO D4 / alone toy 1, 3F20655/910X_,'N-C7 
ee 6 Or oe tes 
L1OX, 'N-C8 Ole 2 Ue oie f 9 Os ate Pe ra) 22 OPPs 
REA. 206) CY MONE (CT \ig tat 56) 
ZOO FOURMAIT SFr 1 5.8 } 
WRT Ge 215) 
ZU ST PORMATICEL 5 SSIS ST Toei ms MOLE ERAT TN Mr/1) 


G COMPOSTTION OF COMPONFNTS MUST BE ENTFERFD IN 

G Toe SAME RORDER AS oe OUm: Pi een RM AT 20 2 

C 
Weer as 207 CY MOA) 6 l= erlor) 

20 7 FORMAT (10X5" COZ Ve 206 Ds / 6 LOXs UNG ie ay, 

xe LOX, H2S ', 
VEZO so. 7 6 LOX > UCL ieee SO es 9 LOM, UCZ 2 Dis / 
ty EO) Xe © C3 ', 


PE2OS Ss /l OX 9) 1-65 (,F20.59/310X5 *N-C4 Ge Oc ice, 
ety Ape G ', 
VePeO. 5.7 >, LOXs "NES e206 5 of eh LOKs 1 =—C6 PE Pie el 
*%~ IDX ,'*N—-C6 ty 
VE 20s Sey LOK LC Vee O sD ey Ws sm Ge 1 P20 co ad 
*,y 1OX_,'*N-C& Us 
WEL ate 7 5 Leg ie A 2d (ele eal Ay Ay aa Aa) 
RFAN(5,206) DIAM,DEPTH,T SOD 
REAN( 5,206) QyTSDyPSD,TDD 
WRAL eV tity 2d, oo 
Do FORMA TC / 7+ LOX, PRODUCT LON PROBLEM 
WRTTE(G, elo) NITAMsDEPTH 
WRITE U6, 207) O,TSOD sTSD, POD. PSD 
DiG FORMAL OLT LL felony LOR UNG I soot gt he olden 9 | 
ee EM CHES 7 5 


110%, "THE LENGTH OF FLOW STRING Sie e rons a 
FFET!) 
217 FORMAT(10X,"FLOW RATE ',30X_'="4F10.3," MSCF/DAY's/y 
110X,'THE AVE, TEMP. MF THE GROUND SURFACE ='yF10.3,! 
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MAIN PROGRAM sia (COM YD) 
LOY aoe TEMP, OF FLOWING STREAM AT SURFACE =',F10.3,' 
rapa TEMP. OF THE PRODUCING FORMATION =",F10.3,! 
i] 
ee WELL HEAD TURING PRESSURE Sie ee Oe aig 


PSL AENL7 A) 


WRITE OUTPUT HEADING 


WRITE (65111) 

FOURMAR (LG, // 74 20Xs SENS TT IVITYOANADY STS OR PAVE T AND 
LyMEBHON t) 

WREkEw C6511 2) 

POURMAT (/7% 18X93 “DEPTH, 12k "0's exe P 25 LSK 5 YT 1 3x 


Ky" T2",4X, 


TBO iIOM HOLE PRESSURES 2X CHANGE. 7.7) 
(a= OPT tA 

X( 2) =O 

X(3)=PSND 

X(4)=TDD 

X(5)=TSD 

NARIGN=X(K) 

DOT LO re y22 

AK S16 Tet CS PLA NOG Sh 4O' Re RGM 
ie deeb) 122 oh 
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CAVES SUB RON TI NEP ROPs FOREACH THE 
PSFUDDER ITT CALA TEMPERA TUR BR ANDMPRE'SSURE? OF GAS 


CALL PROPTAVEMWS Z'PRECs77TEMCsYMOLF>YMW, TEMCsRPREC, SG) 
TRESTPM/ZTEMC 
PRC=PRE/ZPREC 


CB ees) BuNO) TDN Zag etl cere econ ill | tee ae tp aT 


CRE ACh TCT RE a PRE Ze AD 
NENS=0.00149256%*AVEMW*PRE/(7* TEM ) 


CAEL SUR RUT. Neen tC OO CAC WA ir CA SV SO Se Ty. 


CALL VISCON(TEM,PRFE sAVEMW,DENS »VISM) 
RFEN=20,.094595%0%*SG/(VISM*DTAM) 


CALL SUBROUTINES ERTL CT rar CACWwR ATES ERI CTIONS FACTOR 


CALL FRICT CES DIAMS Fie Rte ) 

S=7. OX SCRNEPTH/( 53.34%) EMEZ) 

T= (PSDx*2)* (EXP (S)) 

Wea? 5.0% ((07 10005 8) #2 De SGeT EMER DEP TH*Z 
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PELP=ARS CPREI=PREZD 


COMPAR EOCACCULA T EDS PIR Isr sens sie) sad 


TE (NFL P=0.1) 556,556,560 
PRFIL=PRF2 
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ME h= ore Perr ooG 
FRR=((PREI-DIF)/NDIF)* 100. 
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* THIS PROGRAM CALCIILATES THE EFFERCT OF > 
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* = THF CALCULATED BNTTOM-HOLF PRESSURE. * 


x X(K)=1 VARIATION IN DEPTH x 
x X(K)=2 VARTATION IN FLOW RATE xk 
X(K)=3 VARIATION IN WFELILHEAN TUBING PRESSURE * 


* X(K) =4 VARTATION IN BOTTOM-HOLE TEMPFRATURE * 
*K Ri KJ S5eVARTATION IN WELLAEATD TEMPERATURE < 
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NOMFNCLATURE 


VMOUR 3ar7 MOLE CERAGH LON 


Y MW = MOLECULAR WEIGHT 

TEMNE€e lee CRITICAL TEMPFRATUREs DEGREE R 

PREECE t="CR UTICA SPRESSURESPSLA 

DDAM == 2INS IDE DIAMETER ODF FLOW PIPE. INCHES 
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TRC ScRERUCED TEMPFRA TURE 

PRC =(RERWUC EOS PRESS WR 


ZTEMC = PSEVUNOCRITICAL TEMPERATURE,NFGRFE R 


7PREC PSEWDOCRITICAL PRESSWRE,PSIA 
SG = 9GAS GRAVITY 

HENS a. =. DENS TTY 

REN RFYNOLNDS NUMBER 


AVEMW = AVERAGE MOLECULAR WFIGHT 


VESitg =e WS COS LT Y eCP 
TENMG®. =, CRIT TCANTIPEMPERATMURE SDEGREENR 
PREG. F==CRITICAL PRESSURE, PSTA 


NIMFENSTON YMOILF (16) 4YMW(16) 4 TFMC(16),PRFC(16)5A(8) 


Ke X(5) 
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X(5)=TSD 
DEPT K= 5 5 
NARIGN=X(K) 
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AT AS “A Ih ea RF 
1 X(K)=NRIGN 
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N=X(2) 
PSD=X(3) 
TON=X(4) 
TSN=X(5) 
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TAN=TNN+459.67 
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TRIAL AND FRROR FOR BOTTOM HOLE PRESSURE 
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MAIN PROGRAM eas (CONT) 


CALL PROP( AVEMW,ZPREC »ZTEMC,YMONLF sYMW,TEMC,yPREC, SG) 
TRCI=TND/ZTEMC 
BRCM =rRel/ZPREC 


CABER STERN T TNE Z CAT Ame eA CIT Eee ACTOR 


CALE VARATLUURG WRC leyi/alry A) 
PREZ=) SD LieEMc 

PRF2=PSN 

PRCZ=PRPEZ/ ZF REC 

CARER ZKATZ(TRE2,PRC2q7 7 5) 


CAECUM ATE AVERAGE DENS It x 


DENS 1=0.00149256*AVEMW*PREI/(Z1* TDD) 
DENS2=0.00149256*AVEMWPRE2Z/(Z2*TSD ) 
DENS =CIDENS TEDENS 2 )7 220 

TEM=CT DDRTSD)/ 7. 

PRE=(PREI+PREZ) 72.0 


Chi  SURROUT LINER SCOS TOR CAC A Eerie GAS VIS COS Ty 


CALL VISCN(TEM,PRE sAVEMW DENS »V TSM) 
REN=20.094595*0*SG/(VISM*D1TAM) 


CA =SUBRIL UNE SER GT T DIGACHIES TE rR gt DONG ACT GR 


CAME GT “Ces DI AMSIE Ey REN) 
§=0 .0093 74*SG*DEPTH 
TePRED/t (Os 21 PST SheeZ 2 ) 


PRE3=PSD+( S*T )+V/ (UT ) 
DE WP= A RSP RE 3 PRs) 
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TF (DFLP-0.1) 55645565560 

PRFL=PRF3 

GQ TN 555 
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FRR=((PREI-DIF)/NIF)*100. 

WRITE (69113) DEPTH: 9PSN_yTIN, TSN yPRELyFRR 
FORMAT (19X,7F15.5) 
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MAIN PROGRAM 
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* SENSITIVITY ANALYSIS OF CULLENDER AND SMITH * 


* METHON x 
* THUS PROGRAM CAGE ATESMI HE EFFECT OF * 
VARTATIONS IN THF INPUT PARAMETERS ON > 
* THE* CALCULATED BOTT OM=FURE SPRESSOURES = 
— AOR) =) VAR TAT he uN Se Pn * 
x X(K)=2 VARIATION IN FLOW RATE ** 
x X(K)=3 VARIATION IN WELLHEAD TUBING PRESSURE 
me X(K)=4 VARIATION IN BOTTOM=-HNLE TEMPERATURE * 
* ACK) =5°VARTATION IN WELLHEAD TEMPERATURE 8 
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NOMENCLATURE 


YMOLF = MNLE* FRACTION 

Y MW = MOLECULAR WEIGHT 

TMC w= CR) TICALs TEMPERATURE. DEGREE R 
Crete ent CAL PRES oO RE yi he 

PAM = INSIDE DIAMETER OF “FLOW PIPE, 1 NCHES 
jesy?, = WELLHEAD TEMPERATURE SNEGREE™ F 
TON = BOTTOM—HOLE™ TEMPERATURES DEGREE F 
PSD = WELLHEAD PRESSURE ris VA 

Q = FLOW RATF 

TRC =e REDUCED SFEMPERATURE 

PRC REDUCED PRESSURE 


LREMC =£PSEUDNERITICAL TEMPERATURE, DEGREE R 
APRECs = PSEUDOCRITICAL® PRESSURES PSI“ 


SG ="OAS GRAV Dilby 

NENG] f= DENSITY 

REN =" REYNOLDS SNUMBER 

AVEMW = AVERAGE MOLECULAR WFIGHT 

Virsa = VESCOSTLY.Ce 

TEMG 2 =" CRITICNE TEMPERATURES MIEGK EPR 
PREGA = CRITICAL PREeoUE E>? sia 


NDIMFNSTON YMOLF(16),YMW(16),TEMC(16),PREC(16),A(8) 
»X(5) 

EST Aw AY 0631 SOC 2ST eh. O46 70997,0 4078327 2950.93 5301714 
0 .61232032,0.104888) 3,0.68157001 50 .68446549/ 


RFAN AND WRITF INPUT DATA 


WRITE (6,208) 
FORMAT('1') 

REAN( 5,201) (YMW(T),TEMC(I) »PREC(I)yT=1y16) 
FORMAT(3F10.5) 
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MAIN PROGRAM eee (CONTIN) 


21am FORMA S/S f/f 29% 5 §MOl o WT) artOXs CRI Te. TEMPS ie OX 
SouuG Re laitee fRitG Seer tet 5 
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MAIN PROGRAM eee (ONT 4 Dy) 
* This 7s 
LlOKE TRE: TEMP OF THE PRODUCING FORMATION =t, F110 .354 
Pd Bay 
L1OXKe THbe Whi HEADe TURING PRESSURE IA re Ohne 


Soe ae Ais / 7°) 


WRITE NUTPUT HEADING 


WRITE (69111) 
111 FORMAT ('1',///,30X,'SENSITIVITY ANALYSIS OF CULLENDER 
*AND SMITH MEYHOD # ) 
WRITE (6,112) 
LY 2 GORMA Tay le/i/ey V8 Xie OE Piet yale Xe BOLUS eX at D2) 1 aK syle, 118% 
mia, UA 2 a 
1'BOTTOM HOLE PRESSURF',2X,!' CHANGE ',//) 
Xe) DE PITH 
X(2)=0 
X(3) =PSD 
X(4)=TNDD 
MOSH TS 
Om IKE 1 
NRIGN=X(K) 
MONLO Kia, 92 2 
X(K) =(1.+( J-11)/100.)*NRIGN 
ape 2 251 
1 X(K)=ORIGN 
GOP Tel eal). 
2 NEPTH=X(1) 
N=X(2) 
PSN=X(3) 
TDD=X(4) 
TSD=X(5) 


BORGOMRHOLE PRESSURE  CALGULATION USING CULEENDER AND 
SMITH MFTHOND 
Si NGA IWOeS TEP CALCUL AT LON (SCHEME 


GN1=0.9 
NEPTH=0.0 
TEM1=TSN+459.67 
PRF1=PSN 
SG=0,. 64 


CALL SUBROUTINE PRORSTOSCACUL ATES Tri 
PSEUDOGRIT ICAL. TEMPERATURE AND PRESSURE DF GAS 


CALL PROP(AVEMW,ZPREC »7TEMC sYMOLF »YMW,TFEMC,PREC,SG) 
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MAIN PROGRAM oe J UCONT YD) 


TRC=TEML/ZTEMC 
PROS=PRELAZPREC 


CAS SUBROUMINE ZKATZ. TO GACULATE Z=FACTOR 


GAAUBITPZIKAT ZC TRC e PRC» 25 A) 
N=PREL/(Z*TEML) 

R=NDex? 0 

TEM=TRMT 

PRF=PRF]1 
NENS=0.00149256%*AVEMW*PRE /(7%*TEM ) 


PAIS BROT ERS PO CACULATS THE GAS VISCOUS LY 


CALL VISCO(TEM,PRE,AVEMW,DENS »,VISM) 
REN=20.094595x*0*SG/(VISM*DTAM) 
F=) .0006 


CAITUYSUBRDUTIN ATARI Git 2 C AGWOA Te PARI CTTON FACTOR 


CAI eERTCT (EE sDIAMS Fr. REN) 

Ca(66586" ((07 1000500 42.0) FF ly (bP AMasS. 0’) 
GN=N/(B+C) 

GN2=GN 


GhiecumAr e 1 AT WELLHEAD CONMLT TONS 


M=?. 0*SG*NFEPTH/53 .34 
W=GN+GN1 

NFLP=M/W 
PREN=PSN+NFLP 
GN1=GN? 

hE Pala PU are 
M=2,0%*SG*XNFPTH/53 234 
W=GN+GN1 
DFELP=M/W 
PRFEN=PSN+NFELP 
NTFF=ARS(PREN=PRE1) 
TF (NIFF-0.1) 560,560,555 


CALCULATE IT AT HALF=DFPTH 
PRF1=PREN 


TEM1=459. 67+ (TONNF+TSD)/2. 
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MATIN PROGRAM oferet (CONT ND) 


PRC=PRF1/7PREC 

TRE=SHEM LT AZT EMG 

CALL 7KATZ(TRCyPRC_Z,A) 
N=PRE1/(27*TEM1) 

RaNkx*2?.0 

PRES GPRE LPS) )/ 27.0 

TEM= TFM1 

NENS=0.00149256%* AVEMW*PRE/(7*TEM) 
CALL VISCO(TEM,PRE ,AVEMW,NDENS,VI SM) 
REN=20.094595*0*SG/(VISM*DTAM) 
CAVESERICT (E, DIAM, RE REN) 

C= (66606 ("(0/1000.0) 42.0) FFF) / (Ot AM 5.0) 
GN=N/(R+C) 

M=2.0*SG*NFPTH/53.34 

W=GN+GN2 

NFL P=M/W 

PRFN=PSN+NFLP 

NIFF=ABS(PREN=PRFE1 ) 

PROD UR RO) EP a60re 5605 555 


CARCULATE TAT BOT TOM=Hile CONDE Te hNes 


560) NEP THSk( 1) 
PRF3=PREN 
GN3=GN 
M=2. 0* SG*XNEPTH/ (2 .0%*53.34) 
W=GN+ GN3 
NFLP=M/W 
PRFN=PRE34+DELP 
NIFF=ABS(PREN-PRE1) 
[Ee OLEP—Oek)? 57055 (05354 

554 NMEPTH=xK(1) 
TEM1=TNN+459.67 
PRFI=PREN 
TRC=TEM1/ZTEMC 
PRG=PREL/ZPREG 
CAIIS ZKATZ(IRG,PRCsZ oA) 
(ee / a7ce Te Mis) 
Rae 2.0 
TEM=TFM]1 
PRF=(PRF14+PS)D)/2.0 
NENS=0.00149256%* A VEMWXPRE /(7%*TEM ) 
CALL VISCO(TEM,PRE ,AVEMW,DENS »VISM) 
REN=20.094595*0%*SG/(VITSM*DIAM) 
GALI FRIGT (CEsDTAM,FFsREN) 
C=(666. 6*( (0/1000 .0) **2.0)*FF)/(NIAM**5.0) 
GN=ND/ (B+C) 
M=2. OX SGXENEPTH/ (2 -0%*53 2.34) 
W=GN+GN3 
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MAIN PROGRAM ere © WLAN TT 2at)) 


DEL P=M/W 

PREN=PRF3+DELP 
NIFF=ABS(PREN-PRF1) 

IF (DIFF-O.1) 570 95705554 
CONTINUE 


SE Ss EMERSON” Sta Cate Dri Bo Aa titi F 
BOTTOMS ie PRES SIR 


GN4=GN 

NEPTH=X(1) 

G=GN2+ (4.0*GN3)+6N4 
F=SG*NFEPTH*3 .0*2.0 

P=IP SUF 7/1153 S344G6))) 
OMAP S109 Is LG 
FRIR=(@P=DP LF) /OLF ee roo. 
WRITE (69114) DEPTHsQO,PSN,TND, TSN »P,FRR 
FORMA TH GMOX str 15 25") 
CONT INUE 

CONTINUE 
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MAIN PROGRAM 


ste ste sle slesle ste ste le ole ste ste ste gle ste cle ole ste ste ste ole ste ste gle ale ale gle ste gle sle cle te ale ste ste ste te Sle ale ste ole gle ole ale Sle dle gle sl se ole ale gle 
Sle e es x x Sk deste ste he y 
ASAE TK ASAE FS SS DE SIS TIS OI SE BIS SIE SIC TIS SIR SI FI BS SIS SIK IK SIX SIS DIS SK OSI SK OE SIC IK OK AE SIS SIS SIS SK TK OK OK OK OS OS OK SK KK AS 3S OK 


aN SHEINS: g EVGUTIYS PANIES TS" OOF: SPO TeMANNG «MEAD 6 
“ie PAILS ePROGRAM CAMGUIPAT ES eth naRrRGT « OF a 

VARA LONS IN THe SUN PUTS PARAMETERS (QN a 
* TAP SsCALGUIGA TED BOTTOMS HOE se RESISUR'EZ. * 
x ACK) =1 VARTATION IN DEP TA a 
mK X(K)=2 VARIATION IN FLOW RATE a6 
* MOK} =SReVARTATION IN WELLHEAD TUBING PRESSURE * 


3 X(K)=4 VARIATION IN BOTTOM-HOLE TFMPERATURE * 
3K AGK=S VARIATION IN WELLHEAD TEMPERATURE 7 
Bish Ad ma Paes DAG aP aed bHob le Sand oR nord TE Ad mld pwd pied apa PHo MAb itd oho br aha pa Apyabr ip dpdd be Ad Chg Hee IAEA N bud Ei phd bed ba dou apadpnd ~idpwd ld 


NOMENCLATURE 


Y MOE = MOLE (ARAC THON 

YMW = MOMEGULAR WE TIGHT 

RhEMG Et. CRIBICAL TEMPERATURES DEGREE ik 

PREG AT "OR LVICAY /PRIESSURE SPs PA 

EAMES = SINSL DE DIAMETERS Der EOW PIPE TNCHES 
LSD = WELLHEAD TPEMPERATURES DEGREE NE 

TDD = —BOTROM=HOLE MEMPERATURIE SOEGR BE) 
Psp = WELEHEAD PRESSURE. STA 

-) =e iUW RA lis 

WRC = REDUCED LEMP RR ATURE 

PRC = (ReVUCEI Me RES SURE 

AVEMG *4FePSEUDNDGRITICAL TEMPERATURE sMPEGREE R 
PRRE Cat sPSEUDGCR IT 1CALs FPRIESSURET PS TA 

SG = GAS GRAVITY 

DEIN. @ a=) (OE NSP TY 

REN =—REYNOLDS: NUMBER 

AVEMW = AVERAGE MOLECULAR WEIGHT 

ESM SPS COSRTY »CP 

TEMG. % = IORI CA 7 EMPERIAT URE EG Rise eR 

PR» “eC 1 ICA PRESSURR Pod A 


NDIMENSTON YMOLF(16),YMW(16),TEMC( 16) ,PRFC(16),A(8) 


4% yX (5) 


AT AWAA/ ONS S 0 G23 ORL S4 1099 pO S78 37'7129.0. 53530111, 
1 0.61232032 ,0.10488813,0.68157001,0.68446549/ 


RFAN AND WRITE INPUT DATA 


WRITE(6, 208) 

POR IM Aan LF") 

READ (5,201) (YMW(T),TFMC(1T),PREC(1),1T=1,16) 
FORMAT(3F10.5) 

WRITE(46,212) 
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MAIN PROGRAM eee (CONT'D) 


Ae late PRESS oi! 8et/ 5 
14 7X, "DEG. RANKIN' 412X,_ *PSIA',/) 
WRITE (6,202) (YMW(1) ,TEMC(I) »PREC(1),I=1,16) 


FORMAT (CLOX4 §CO2 ",3F20 6597/5 10X_'N2 de ee) ad Oca ye / 
Ea MO >, QPEL WSS U 
Ve SRS 0. os / 5 LOX. BEE Me TIA) e DIS fas Ne Ae 153F20.5% 


ed gs LA on 

1 PSE 20. Sp f5 LOX sg U1 = C4 1, 3F20.59/510X_'N-C4 ; 

* 9 3F200.59/910X_t*I- 

ee eH aren es) Oy el Ce eae Mic Heo AO. e/a yt 1 Co { 

Mees Dig / 9 LOX 

1LYN=C'6 by DIP oD 9 et es LO, ty OE ZO0 ey / 4 LOK N=CT 

Ar Or 2 Ores 1% 

PIOX, 'N-CS "y3F20¢59/,10X_'HE iy O20.) 
READ 5s 206)" (YMORECT) .l=15 16) 


ORM ATECS FT 5). 3)) 

WRITE (6,215) 

PURMAM OULU SSS LISI TT ved Xe (MORZESERACT CON V5 } 
GOMPOSTTION OF COMPONENTS MUST BE ENTERED IN 

tate oh ROR OER AS THe ent PUN -EORMAT 220.2 
WRITE(6,207) (YMOLF(I),T=1,16) 

FORMATS «( LOX4¢CO?2 "sF20 659/79 10Xy_ INP te Eee Dat 
*%_ LOX,"H2S ty 

DR ZO". 29 LOK SE Cl VF 20 Def 9 LOX CZ Hy F2Z20655/ 
*, 1) X,'C3 'y 


PP20 35405 VOX,' T=C4 "sF20 6559/9 10X_ (N-C4 ie 20 eee, 
Kee Xs OL =—CS ', 

1F20.59/,10X,'*N-C5 F209 LOK, Hl =C6 Wy eZ Oe Diy 7 
Ky 1OX,'N-C6 ', 

VEZO Ds / pvOXs OL = Ce U HIP ZO 6 DAH OOK G5 UINSIG 7 Ore itaiel Oot D) 9h/, 
*~ 1OX,'*N-CB ', 

IP2040457 7 LOK AE Ver ZO ey SSS SS) 

READ (S¢206) OLAM DEPTH FSD 

RFAN( 5,206) QyTSNyPSN,TDN 

WRITE(6,219) 

FORMAT(///4910X_ *PRONDUCTION PRORLEM® ) 

WRT EUs 2 Oye) VAM DEP Tr 

WRITE(6,217) Q,yTSOD,TSN,TOD »PSN 


216. FORMATWAL////slLOXG TURING LT. Doty 2X tat. Flo. 3,4 


He JUNIE SY Ao 


110X%, "THE LENGTH MF FLOW STRING SN ee a 
acu ft Ese) 

FORMAT(LOX,"* FLOW RATE! ,30X,'=',F10.3,! MSCR/DAY * 97% 
110%, THE AVE. TEMP. VPS THE IGROUINOSSUREACE =',F10.3,! 
x Fly /, 


110X,"THE TEMP. OF FLOWING STREAM AT SURFACE ='yF10.3,! 
6 F's, 


Ve ost, 


‘a 
a PAT Peery ” age at 
} . } ai a 
’ ?,/= Pons t ' Pb - 
nau a 
) i }—j eT, .* sO tie yt map. oh 
} . . if z. ; | 
Aw At Ma eh! OF Gs ae ah ei a 


if cea «* ist (iad tdaaaaeene 
bad yo) et aes 
ae 
(\ nae pt hive un Paty OANG sr am 
Be d 


wT +petring an joan 29ers Fe i ut znamern 
GD PAMHA TOO ANY 2h Be yattA2 at) 


. 7 


if 
(al, cede MNAanIME LANDA ADST EAR 
Vee f04 »! LW xT hes BAe! ; ‘ROD? hai 47 sean 
ti ‘KEL *. 
eas weg SIT sAGTS Cet. tol 
«ts 
Ayes inte’ i =4* XOF sXetslhty ! eb teil i , s 
it) fey el) ae 
Vy c0Sy? SOT RM Ret Po mwe His \ e088 
Gg! aan oh OL a 


LdEeUS4s! i.) AY oxinia Wa atts! Toh aie Sue at 


« 


rorgglle a eut 
Sher gpl male ¢ 
Pe ie) 


° 
a) 


amy am, 


Vgk TAY t= ere Not 
id <ney at> : | 


Fe oe if wap 
= _ 7 


tee 


<=) 


>) 


EG) 


Les 


byZ 


88 


Tet 


123 


MAIN PROGRAM eee | CONT 2 DD 
TIO%. Wine TEMP. OF THE PRODUCING FORMAT ION =MgrtO. ae. 
>t Fi, /, 
1tO%. STHE WELL HEAD TUBING PRESSURE SA 1) are" 


Fm Ana, f./ ) 


WRITE. OUTPUT. HEADING 


WRITF (65119) 
PORMAT CML SS fe chs 


*U SENSITIVITY. ANALYSISSOFP,POETIMANN METHOD!) 


WRITE (649112) 
ORM AT Oe 5 LBs UDP RA Lox eo Se gle ky ie 2 od Set teen oA 
get Hie! 5 X's 


L'BATTAM HOLE PRESSURF',2X_" CHANGE',//) 


Kt 1d] Die Pa A 

X(2) =0 

A{3)=PSD 

X(4) =TOD 

X(5)=TSD 

IMMA) Bs K=l,5 
NRIGN=X(K) 

ND 44 J=1,22 

X= Ci.+ (3-11) /vwoo.) OR TGN 
Ce 63= 20 ie Ol (ines 
X(K)=ORIGN 

SOT Oe 33 
DEP Tia = KG ol ) 

M=Xx( 2) 

PSND=X(3) 

TDD= X(4) 

TSD=X (5) 


BOTIOMNMHOGLE PRESSURE CALCULATION USING 
POFTTMANN METHOD 


F=0 .O0N06 

TDD= THNH+459 06 7 
TSH=TSI4+459..67 
TEM= (4 ODT SD) 7 Si 
PRF1=PSD 


CAL SSWAROUT ING TPROP Sh CACULRATE Thr 
PSEWDOGRIDIGALeTEMPERATURE AND PRESSURE DOF GAS 


CALL PROP (AVEMW,ZPREC y7TEMC yYMNLFyYMW,TFMCyPREC,SG) 
TRGE2=) eiAAeAREMG 
PREZ=PRELAZPREC 
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MAIN PROGRAM spate UIN dae wn) 


H=0 .05 

Y1=PR2/H 
N=(2*Y1)/2 
NNIM=N—3 

NO LO 1, NN BM 
PR( I )=0. 2+H* (I-1) 


Chis! MSUBRIU PINE SKA, Bi eOAGIEA TE ZAG TOR 


GALL ZKEATZECTR2,PR CPW, 701) 5A) 
VERS rere PR Ch) 


CAEL *SHB ROUT INESOSE TO EVALIATE THE INTEGRAL 
AE 27 PR NUMERTEALIEY 

be ioral aed Sala Oa ek OO Ns a de Sy en BRINE alice 
NUMFRTICAL INTEGRATION 


CALL QSF (HyVsAIN,NDIM) 
FRI1=PR2-PR(NDIM) 


ASSIGN AN EREGR TOLERANCE (Or On.001 SIN PR 


Le CARS CERT=OCO00W |) 22522921 10 
C2=AIN(NNDIM) 
Gu "55 


WSESTRAEPEZOIDALY RULE TOCUBTAIN THES VALUE OF 
THE INTEGRAL™AT® THE REGULREN PR 


V1=V(NDIM) 

AL=AIN(ND IM) 

CALE PYKATZ EUIR29PR2¢722 4A) 
V2=72/PR2 
C2=Al+(V1i+V2)*ERI1/2. 


TRIAL WAND FRROR SPOR THe BOT OMFG Be PRESSURE 
GUESS AN INITIAL VALUE FOURS TAR Roi OM-SOLESPRESSURE 


PZ= 3 00. 0) 

AP=(P2PPRETY/ 2% 

NT=TEM 

TRASAT/Z1 EMC 

PRA=AP/ZPREC 

CALL ZKATZ (TRAgPRAg/7AyA) 
NENS=0.00149256%* AVEMW*AP /(ZA*AT) 
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MAIN PROGRAM siete (Catal 2210) 


CALS SSUBROUTLINERV DS CO eLOSCAGCUIEAT ATH ERGAS AES GO Silvey 


CALL VISCO (AT,AP,AVEMW,DENS,VISM) 
REN=20.094595%*0Q*SG/ (VISM*DIAM) 


GALE SUBROULINGRERIC TINO ACACUMATEIERICTION FACTOR 


CALL FRICT(E,DIAM,FF,REN) 


W=NT AM%% 5 

PR1=P2/ZPREC 

TR1=TR2 

nS (OS 

Y2=PR1/H 

N= (2% Y2)/2 

NDIM=N=—3 

MO 21 JL=leNDIM 

PR( J )=0. 2+H*(T—-1) 

CARL ZKA RE v UPR P RIGGS Zale), A) 
V(I)=Z(7T)/7PR(1) 

GAEL OS.F SCH VAT Nee NID) 
FR2=PR1-PR(NDIM) 

VR OAB SHER 2=—O1s 0 OR)? 5555 5766 
C1l=A IN(NDIM) 

GO TO 94 

V1=V(NNDIM) 

A2=AIN(NNDIM) 

CoAd lei. KoA TeZer (RRs PAR seh FAY) 
W274 SRR} 
C1l=A2+(V14+V2)*ER2/2. 
CONTINUE 

MS=( 53, 34*A1/SC)2 (Cil-C7) 
TaWwxe( XS%** 2?) JU 

NP=( (NFEPTH*T) /(XS-DFPTH) )**0.5 
PN=PSN+NP 

NFLP=ABS(PN-P2) 

TF (NFLP-0.1) 98,98,99 
P2=PN 

GOTO 296 

h1P= 3635.29 
FRR=((PN-NDIF)/NIF)*100. 
WRITE (65113) DEPTH,O,PSND,TOIN, TSN yPN,FERR 
POR MAS £(HOXs 761555) 
CONTINUE 

CONTINUE 
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MAIN PROGRAM 


He ste slo Sle gle ale Ale Se sie 5! te sf. te al, 4, D A 
l- ste ste le yle cle gle slo gtr ste ste gle ste gle Se cle gle ste ole slo gle she gle gle yey te ste ste ste ste ale slo Sle ste ste Se Sle ole gle cle gle ste Sle gle ale ote 
DE BS BS AS ASAT YS DS BSL DYKE DS IE ALK DK DIS DIK AK DIS DIS DIL DK DK DIE DIS DIK DIE DIK DIK DIS DIK DKK DIK DIK OK BK DR IK OK AK OK OK OK AS 38 OS FS 3K AK OS 


PSENSHELLV INE ANALYSIS Oe SUKKAR TAND CORNER METHOD> 


THIS PROGRAM CALCULATES THE EFFECT OF x 
VARIATIONS IN THF INPUT PARAMETERS ON a 
aK THE CALCULATED BOT TOM=HOLE® PRESSURE. es 


a MEK = VARTATION IN DEPTH ** 
x MK) =2% VARIATIONS IN FLOWRATE * 
x X(K)=3 VARIATION IN WELLHEAD TUBING PRESSURF * 
* X(K)=4 VARIATION IN BOTTOM—-HOLE TFMPERATURE * 
aS OUR) FB VAR TATION TN WELCREAD TEMPERATURE a 
Pen P hd >kS Cierra ph ab Agr tae Hah ha Agama diab idher Oridaparnaenamiapha mime nd niqnraetdphapnae ap htoi ar nariqirlae shit taba redpadpadns dn nara nud cade ara gerd 


NOMFNCLATURE 


YMOUF: = MOLE iARAGTIDON 

YMW = MOLECULAR WEIGHT 

LeMG ee, CRIS CAL TEMPFRATURE, DEGREE F 

PRiCes = CRIT ACALy PRESSUREGIPIS TA 

PLAM == INSIDE DIAMETER Ue FLOW PIP EST NCHES 
YSp =) WrisL HEAD SREMP ERATURE + DEGRIE Gast 

TDD = 8 Vi OM= HOLL Pare MPRA TU Ri DE GR Ee ie 
PSO = WELLHEAD PRES SURE.P S14 

a) = FLOW RATE 

TRC = REDIC Bbw eMPERAMURE 

PRC = IRIFDUGRDEPRESSURE 

7TEMC. B&B PSEUDOCRITICAL TEMPERATURE, DEGREE R 
EPIRECS =? PSEUDOCR IT TL CAIQIPIRESSURET Piss th 

SG = GAS GRAVITY 

DEMS. S.=7 DIES DT ¥: 

REN = RFYNOLDS NUMBER 

AVEMW = AVERAGE MOLECULAR WEIGHT 

VIS ye =eWIS COSTE Y—CP 

TEMG, S27 CRDTICAL TEMPERATURE. DEGREE os 
PRECHS=eCRATICAL PRESSURE SPSITA 


NDIMENSTON YMOLF (16) YMW(16),TEMC(16) »PRFC(16),A(8) 


*_ X(5) 


NTT N7106341 5062370467099 JOS 78392729 »50-.5393077T1, 
1 0.61232032,0.10488813,0.68157001 ,0.68446549/ 


READ AND WRITE INPUT DATA 


WRITE (6,208) 

FORMAT(! 1") 

REAN( 55201) (YMW(1),TEMC(1) »,PREC(I ),1=1,16) 
PORMAT(3F10.5) 

WRITE(6,212) 

FORMAT(/S////429X%o MOL e Whe ty 1LOXy* CRIT. TEMP. '48X 
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MAIN PROGRAM sere cM, ante) 


geet emaie ee S Stents Zils 
LSI XWU DEG 6 (RANKING. 12x PS TAS J) 
WRITE(6,202) (YMW(1),TFEMC(I),PREC(I),1=1,16) 


232 FORMAT @10X ¢'CO2 Veo ee Dal VOKs NZ eB 20s e 7 
Tell) Kg ies : 
Veo Aes Ss VOX MO) Hy 3F 2025974 VOX ' C2 Mr yratac las 


*/ 5 VOIX, * C3 
1 Vy oF 20. 54/5 1LOXs t 1-64 patie) 6 D7 ee N= Oe : 
re Ie A Oe os fe LOX git 
Le Bot 20 eo 9/7. LUA NAC 5 Hy Or 206 Dig fo LOK LOG ' 
ete es fs LOX 
1'N=-C6 UY rot 2 Ole te) gan ON ye ial Cat, Ng Sag 9:/ anol ONen N= Cul, 
x DES el SAO llr ad Lor 
110X,'*N-C8 Ne SCO eee wa kas eee sD LO so) 
READ CS. ZUG) (YMOER CIs 1 = fsilG) 
296 “FORMAT CSF 15.3 
WET EUG 2) 5) 
21D FORMAT CUVE LA, LATS pet kay eer RAG TL CNet s2/*) 


C COMPOSITION OF COMPONFNTS MUST BF ENTFRED IN 

C THF SAMF ORDER AS THF NUTPUT FORMAT 202 

& 

WRITE (6,207) (YMOLF(I),T=1,16) 
207 FORMAT, (10X',* C02 ae Oise) tas IOP iN i P20. oy 7, 

%~ 10 X, SH2S Vy 
VECO el OK scl LOO sao Ge 1,F20.5,/ 
ESO, 8 CS ', 


TP20 5s 7/5 LOX,s' 1-C4 1 ,F20 059 /,10K_'N-C4 ire FeO oioen/ 
*,10X, '1-C5 Ue 
1F20 .59/,10X%,'N-C5 VR ZO eo eels be Gs Wee Oe Oe, 
*,10X, *'N-C6 ies 
LEO. Desig LOX s! 1=CT. SO eos OMe peep yey Ae pO) car pel Fe Oy f te elo 
*%, 10 X, *'N-C8 ', 
LEZ son 7 LAs Ee Dei es OSS ey A ar TAT Ae 
REAN(S, 206) DIAM,DEPTH,T SOD 
REANUS. 206). QO.) SO gh SIs 1.00 
Wie de i eaao see clue’ 
219 ENRMAT(/// ,10X, PRODUCTION PROBLEM! ) 
WRITFE(6,216) DIAM,DEPTH 
WRITE(46,217) O,TSOD,TSN,TNN,PSD 
F116 CORMATIS///7// 5 10%e* TURING In De',27Xet=',FlO.3,! 
INGHES! 57: 


110X,'THF LENGTH OF FLOW STRING =!',F10.3,! 
RF FETS) 

O17 FORMAT(LOX,"FLOW RATE", 30X,"='5F10.3," MSCF/NAY's/, 
1190X,'THE AVE. TEMP. OF THE GROUND SURFACE =',F10.39! 
xe eal ee OF 


110X,'THE TEMP. OF FLOWING STREAM AT SURFACE ='yF10.3,! 
a Fl, /s 


ee OTE! 


| — 


sean den Rem 


AeA ROLY Moe eee tt | ne 
eal i 0 eh f Les ; 
Ye eG Ay] 
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MAIN PROGRAM eee | GON TD) 
PO Ay Ue LEMP... OF) THE PRODUCING «FORMAT LON =',F10.3,! 
Fly/y 


PIO As (AWE OWELL HEAD TUBING PRESSURE 
BS /'/) 


WRITE OUTPUT HEADING 


WRITE (6,119) 
PEO OR Mae (fey e710) Xe 


SENS ERIVITY ANALYSA Sit SUKKAR-CORNE UE METHOD ) 


WR Te 0's. Lids) 


TS 2e EUR MAR O/ /ue LS Xe UP Tita egal 2 Xe UG) a 2 Ke IP Pg WK gallate aK 


PUN Bieae Bey. OF 


L'BOTTOM HOLE PRESSURF',2Xy! CHANGE',//) 


Xl), =DE PLA 
X(2)=0 
X( 3) =PSD 
X(4)= TDD 
X(5)=TSD 
NY: 33 Kals5 
NRIGN=X(K) 
NO 44 J=1,22 
RSC (J) OO) ae LGN 
lees Gelert eas sales 
88 X(K)=ONRIGN 
Cle EO 33 
tt SEPT hx 1) 
M=X(2) 
PSN=X(3) 
TONN=X (4) 
TSN=X(5) 


ROTTOM-HNI_E PRESSURE CALCULATION USING 
SUKKAR AND CORNELL MFTHOD 


TNN=TNN+459.67 
TSN=TSN+459 267 
TEM=(TDD+TSD) /26 
PRFI=PSN 

AT=TEM 


EXCL SUBROUTINE PROP RTO 2G ACULA ICs SH 


PSFUNDNCRITICAL TFMPERATURE AND PRESSURE 


OF GAS 


CALL PROP( AVEMW,ZPREC,ZTEMC yYMOLF »yYMWyTEMC PREC, SG) 


TG eine 7. Altes 
PRC ee i 1, bie G 


COORHT ae sy e 
NEE MET TYNES ys ene THAN 9 
(ber sabato” As 
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tik 


Que GG 
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MAIN PROGRAM ere o UCONN T thy) 


FF=0.01482 
B= (666.6% FF*((Q0/1000.)**2)*(AT#**2) )/( (DIAM**5) *(ZPREC* 


sei ste: ? ) ) 


H=0 .05 

Valea Ree 
N=(2*Y1)/2 

NA TM=N=3 

NO 10 JF=1,NDIM 
PR(T)=0. 2+H* (I-1) 


CELE SUBROUTINE ARAT 2 TO sCACULS TE ZF ac ror 


CALE ULKAL OPTRA SPR (It )is 77 C8 hi A) 
TECIO SA CTP RC ye Croce ea udesem iG Lal) seers, 3) 


CARE SUBROUTINE GSP et *EVACUATE (THE SIN TEGRAE 
OP CL/ PRY (let (BS CZ7 PR 2) NUMER DCA LILY, 
Decor oo Cle le OOo nh PRE eo aust PN 
NUMFRICAL INTEGRATION 


CALL OSF (HeV,AIN,NDIM) 
FRIT=PR2—PR(NDIM) 


ASS JGNEAN ERROR: TOLERANCE THEY O.O0G1) TIN EPR 


PPARs Or R LO Olay) 22 Zin oredr. 
C2=AIN(NDIM) 
Ci LoCo 5 


USE TRAPEZOTDAL RULE TOSURTAIN THE: VAEUE OF 
THe tN LEORAL AV tnt REGU ERED PR 


V1l=V(NDIM) 

AL=AIN(NNDIM) 

CMI EZRATZ (TR2¢PR2Z 9229.4) 
Wi2= E27 PR Z 
C2=A1+(V1+V2)*ER1/2. 

We (SG*NFEPTH)/(53-34*AT ) 
AINL=C2+W 

TRelst eM, 2 hesiG 


AN FRROR TOLERANCE OF 0.00001 IS ASSIGNED 
TN MATCHING THE CALCULATED INTEGRAL 


F=) OOOO] 


ver aaa-5 ied 


| et 
seston elle 


Ane ITM AKT 4 aca OF a ‘o 
IIA ITAAMI Fad BSP 
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MAIN PROGRAM eee CONT UI) 


ND TM=NDIM+4NDIM 

NO Ll T=1,NDIM 

PRACT =e 2a (al 1) 

CAIES ZICATIZIN( TRL PRMD tZ 1% AD 

VICE) =A TSP RICEO) P00 Se (Be ZT RP RAL) ) ea 8) 
CALL QSF (HyVsAIN,NDIM) 

GOATO = 

NDIM=NNDIM—1 


USE He METHOD. OF SENTRA FAG Vy ING STOUR TATE Tia 
VARWERD RT HEMRRE THAT. GIVES@ ANTSINTEGRALPWALUE 
FOUA ES TT) THE CALCULATED INTEGRAL VALUE AT THE 
AO TOM-HOLE COND TT PONS 


LEY ( AIN( NPHM) —Ad NIP ADs ly 2 
GOrqtO. 3 

XH=A IN(NDIM+F1 ) 
PH=PR(NDIM+1) 
XL=AIN(NDIM) 

PL=PR(NDIM) 

AM={ eX KEY) 726 
PM=(PH+PL)/2. 

TF (ABS CUXM=—AINI)—E) 595,14 
TF (XM-AINI) 945,6 

XH=XM 

PH=PM 

GOTPOC20 

XL=XM 

PI=PM 

(Sa) WG) 240) 

P=PM 

PR=P%x 7 PREC 

ORRSLOTaS4 129 
FRR=((PB-NIF)/DIF)*100. 
WRITE (69113) DEPTH,O,PSN,TNN, TSN,PR,FRR 
FORMAT (1 OX, +ESS5) 

CONT INUEF 

CONTINUE 
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SUBROUTINE PROP 


SUBROUTINE PROP(AVEMW,ZPREC,ZTEMCyYMOLF 5 YMWy TEMC y PREC 


* 9 SG) 


Tris SoVBROOT IME (CAR CUEALES SHE PSEURU CRT LEGAL 
TEMPERATURE AND PRESSURE OF THE GAS MIXTURE. 
AZIZ AND WICHERT CORRFCTION TO THE STANDING-—- 
KAT? “CHAR TE TSP USENGROPGRRRECG ENRIMIHE PRESENCE 
OFM SBURE GASES 


NIMFNSTON YMOLF( 16) 5YMW(16),TEMC(16),PREC(16) 
M=]1 


AK SON HE A) eee Ome) aed 

PEeM=0, PSEUDOCRITICAL PROPERTIES ARE CALCULATED 
BY Tt SPECIFIC. GRAVITY CORRELATION BY BROWN AND 
PIAL. oHURT GASs 1S, CORRECTED BY GARRES METHOD 
ANDEWTCHE RT?! ANDeAZTZS CORRECTION FACTOR 


Pee Means PSEUDDCRITICAL PROPERTIES ARE VGALGULATED 
USN GY MOLAL. AVERAGES AND CORRECMED POR SGUR GASES 
WSING “AZT 2 CAND  WECHERT CORRECTION 19 STANDING- 
KAZ. CHART 


TF (Ml) 65494 

AVEMW=0.0 

PTFMC=0.0 

PPRE C=0.0 

MOL Ss tla k6 

AVFMW=AVEMW+ YMOLF (1 )*YMW(T ) 

BT EMC =PPEMCHYMOCE (CI Y*TEME (1 } 
PPRFC=PPREC+YMOLF (1) *PREC(T) 

SG=AVFMW/ 28.97 

COL STi tS 

Pate MGs Wie liaise fae 5G 

PTEMC=PTEMC=8 0.0 SV MNLF (1)—=250. *YMOLE(2)+130. *YMOLF(3) 
LFeriSG=0.85) 7e798 

PPREC=695.1-40.*SG 

Gat 70) LOS 

PPREC=704.396-51 4 724*SG 
PPRFC=PPRFC+440.*YMOLF(1)-170.¥*YMOLF(2)+600.*YMOLF (3) 
SL=yYMOCE CIOS YMOLR (3S) 

S2= VMOU (sa) 

F3=1 20. 0%(S1**0.90-S1**1.6)4+15.0% (S2**0. 50-S2%*4 ) 
S4=PUEMCES 2 O52) ro 
7 EMC= PPP HM Greg 

7PREC=PPRECXZTEMC/S4 

RETURN 

FND 
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SUBROUTINES ZA te 


SUBROUTINE ZKATZ(TR»PR,Z,A) 


THIS SUBROUTINE CALCULATES THE Z-FACTOR USING A 
RENDUCFD B-W-R EQUATION. THE B-W-R CNEFFICIENTS 
WERF CALCULATED FROM THE TABULAR DATA IN 
"HANDBOOK OF NATURAL GAS ENGINFFRING!. 


NIMFNSTON A(8) 

Le eek 

PRC=PR 

a Soa g) 

Th oO) 

IPCTRE=). 05-43 135473 

TP(TRO=3.0) 14,124,438 

TECeRG=—30'.0) 15.15,43 

HOF 2 OO PTER=f, 10 

NR 2=DRx*DR 
T1l=(A(1)*TRC-A(2)-A(3) /TRC%* 2) *DR 
T2=(A(4)* TRC-A(5)) *®DR? 
T3B=A(5)*A( 6) *NR**5 

T4=A(7)*DR2/TRCO**2 

T5=A(8)*DR2 

T6=FX PA-T5) 
P=(TRC+T1+ T2+ 73) *DR+T4 *DR*(1204+T5)%*T6 
NP=TRC+2.0* T1+ 3 OX T24+6 OX TA 4+T4*TOX (3.04+3.0%7T5-2.0%*T5 


NR1=DR-(P-0 .270*PRC) /NP 
PEM RIy 2167.64 le 

DR 1=0.5*NR 

TECDORI=2.2) 19,19,18 
DR1T=NR+0.9% (222-NR) 
TF(ABS(NR=DR1)-—O200001) 21,70,20 
NR=NRI 

Lae 2h PRG/ (PRT RO ) 
CONT INUF 

RETURN 

FND 
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SUBROUTINES NeESCn 


SUBROUTINE VISCO(TEM,PRE s,AVEMW,sDENS,VISM) 


THES SUBROUTINE GALCULANES® THE: VISCOSITY pOR 
THE SGASSMIXTURE 


VK=(9.4+0 .02% AVEMW) *TEM%%*1.5/( 20920419 .0O*XAVEMW+TEM) 
VX=3-5+986 .0/TEM+0.01* AVEMW 

VY=?.4-0. 2% VX 

VY 1=VX*NFENS#* VY 

VISM=VK*FXP(VY1)/10000.0 

RETURN 

FND 
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SUBROUTING PRICT 


SUBROUTINE FRICT(E,DIAM,FFyREN) 


THis ssUBROU TING CALCULATES. THe eR TOI OND rac Tur 


NIMFNSION G(100) 
C=ALO0G(10.60) 

RR=F/DIAM 

VEC Rees oes el 
RENCR=3200.0/RR 
TF(REN-RENCR) 28,32,327 

= Ore Onis 

MSO = 14-1 O00 

AVERR7 Sie (2651/7 (REN*X) 

MZ A LOCK XU) 2—302585 
FX=1.0+2.0%X*X2 
FXN=2.0%X2—5 .02/ (REN*X *X 1*C ) 
G(1)=X-FX/FXD 

X=G(1) 

TFC I=T) 30530,29 
TPCABS(G(T)=G(1I=F) }=0 2000007) 312,31, 30 
CONTINUE 

FE P= Kor xX 

GO-TO B3 

PE=1.0/ (02. 0* ALNG(3..77RR)/2. 302585 ) ** 2 
RETURN 

END 
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SUBROOUTINGS OSA 


SUBROUTINE OSF (HyY¥9Z sNDIM) 


THIS SUBROUTINE DGES “THE NUMERT CAL INTEGRAT LON 
OF ANGIVEN “BUNECT ION WS ING SIMPSON*S RULE AND 
NEW LOIN Ss 9376 RUCE 


DIMENS TON ¥C6)%52. (5) 
HT=.3 333333%*H 
TF(NDIM—5) 7,8 51 


Te NDIM SIS GREATER nam 6 
PREPERATIONS OF INTEGRAT LON S200 


SUM1=Y(2)+Y(2) 

SUM1=SUM1+SUM1 

SUM ES et YC SUM + ¥5(63:))) 
AUX1=Y(4)+Y(4) 

AUXL=AUX1+AUXI 

AUXI= SUMIF+HT#(Y (3) +AUXI+Y(5)) 


BUX Dot eA Yl) toe BA Oe WY) CEN AOU) tee aie Br Yet 3) ee 5) 


bic al (oi) ba) 

SUMZ=VICS)+ ¥ (5) 

SUM2= SUM24+SUM2 

SUM? =AUX2—HT* (Y(4)+SUM2+Y (6) ) 
7(1)=0. 

WOKEN (3) + ¥ (23) 

AUX=AUX+ AUX 

TU2ZY =SIM2— HT (CY (20 +AU XY (4))) 
FAS = SIM 1 

7(4)=SUM2 

TFE(NDIM—6)5 4592 


INTFGRATTON LOOP 


Mi 4 Ta 7 gNITEM 92 

SUMIL=AUXKI 

SUM2 =AUX2 

FE) SG i) 1G 2 10 Bes 
AUXL=AUX1+AUX1 

AUX1= SUMLFHT*(Y (1-2) +AUXK1+Y(1)) 
F( 1-2) =SUM1 

TE ( T-NDIM)3 9696 

AU eee) Nr) 

AUX 2=AUX2+ AUK] 
AUX2=SUM2+HT* (Y (T-1L)+AUX24+Y (141) ) 
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SUBROUTINE OSF Sc SSIS EL Be) 


7(I-1)=SUM2 
ZUNE My =A XY 
7(NDYIM) =AUX2 
RETURN 


7 (NOD IM=1)=SUM2 
2 UND My = AUX 
RETURN 


END OF INTEGRATE VON: ame 


DE CNDS Me 4 Velie cliles 8 


NDIM TS "FOUAL TO) 4 OR 45 


SUM2=1.125%HTX(Y (1) 4¥ (2) 4¥(2)4Y¥(02)4+Y (3) 4Y(3) 403) 


a Fee et) &) 


SUMTEYC2Z2)+4Y(2)) 

SUM1=SUM14+SUM1 
SUMLSHT*(Y (1) +SUM1+Y (3) ) 

7s a) mal 

AUX1L=Y(3)4+Y(3) 

AUX 1=AUX1+AUXI 

7 A(Z2) =SVMA=AT = CY (2) FAIIX14Y(4)) 
TE(NDTM=5)10,9,9 
AUDLEY U4) 49 (4) 

AUX1L=AUX1+AUX1 
7(5)=SUML+HT% (Y (3) +AUX14+Y (5) ) 
7(3)=SUM1 

7 (4) =SUM?2 

RETURN 


NODMe iS. OU ARS 1 33. 


SUM DST ret aA LEY C2) vt 2) 5 oe eS ye) 
SUM2Z=Y(2)+Y(2) 

SUM? =SUM2+ SUM2 

PAeo rl ely (ele) SUM2 Ne CB) 

INE Oe 

7C2)=SUM) 

RETURN 

FNTI) 
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